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Abstract  This  paper  surveys  current  approaches  used  in  the  United  States  for 
assessing  casualties  in  simulated  tactical  engagements  between  general-pur¬ 
pose  military  forces  in  conventional  air-ground  operations.  It  first  discus¬ 
ses  the  various  modelling  alternatives  available  to  the  military  OP  worker 
and  then  expounds  upon  both  detailed  Lanchester-type  models  of  attrition  in 
tactical  engagements  and  also  aggregated-force  models  based  on  index  numbers 
(e.g.  firepower  scores).  Methodological  aspects  are  emphasized.  Simple  auxi¬ 
liary  models  are  used  to  illustrate  modelling  points  for  developing  and  under¬ 
standing  complex  operational  models,  but  examples  of  current  operational  mod¬ 
els  that  use  these  two  theoretical  approaches  of  casualty  assessment  are  giv¬ 
en.  Concerning  detailed  Lanchester-type  models  of  attrition  in  tactical  en¬ 
gagements,  simple  auxiliary  models  are  used  to  illustrate  modelling  concepts 
and  issues  such  as  (1)  various  functional  forms  for  attrition  rates,  (Z) 
determining  numerical  values  for  attrition-rate  coefficients,  (3)  various 
operational  factors  to  be  considered  in  attrition  models.  Index -number  methods 
for  aggregating  military  capabilities  and  aggregated-force  models  of  attrition 
are  then  discussed,  and  a  Lanchester-type  aggregated-force  attrition  model  is 
developed. 
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Obersicht  Die  Arbeit  gibt  einen  Oberblick  Uber  die  derzeit  in  den  Vereinio- 
ten  Staaten  gebrauchl ichen  Ansatze  zur  Verlustermittlung  in  der  Simulation 
taktischer  Luft-Landkriegsgefechte  und  -operationen.  Zunachst  werden  die  ver- 
schiedenen  Alternativen  der  Model Iformul ierung  angesprochen,  die  dem  Bearbei- 
ter  mil itarischer  OR-Probleme  zur  Verfugung  stehen,  und  dann  werden  eingehend 
erbrtert  sowohl  ausfuhrliche  Lanchester ' sche  f'odelle  von  Abnutzungsprozessen 
in  taktischen  Gefechten  als  auch  f'odelle  auf  Gesamtstreitkrafte-Ebene,  denen 
Indexzahlen  (z.B.  Kampfkraftzif fern)  zugruridel  iegen .  Dabei  werden  methodi sche 
Aspekte  besonders  betont.  Es  wird  von  einfachen  Hilfsmodel len  Gebrauch  ge- 
macht,  urn  wesentliche  Punkte  der  Model  1 formul ierung  beim  Entwickeln  und  Ver- 
stehen  komplexer  Planspiele  zu  veranschaul ichen;  es  werden  aber  auch  Beispiele 
gebrauchl icher  Planspiele  gegeben,  die  diese  zwei  theoret i schen  7uqange  zur 
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Verlustbewertung  benutzen.  In  Bezug  auf  ausfiihrliche  Lanchester-Model le  von 
Abnutzungsprozessen  in  taktischen  Gefechten  werden  einfache  Hilfsmodelle  ver- 
wendet,  urn  Konzepte  des  Model lansatzes  und  solche  Punkte  zu  verdeutl ichen  wie 
(1)  verschiedenartige  Funktionsgestal t  der  Abnutzungsraten,  (2)  Bestimmung 
numerischer  Werte  fur  die  Abnutzungskoeff izienten  und  (3)  Beriicksichtigung 
verschiedener  operativer  Faktoren  im  Abnutzungsmodel 1 .  Schlieftlich  wird  auf 
Indexzahlmethoden  zur  Zusammenfassung  militiirischer  Leistungsfahigkeiten  und 
Abnutzungsmodel le  auf  Gesamtstreitkrafteebene  eingegangen,  und  zuletzt  wird 
ein  Lanchester' sches  Abnutzungsmodel 1  auf  Gesamtstreitkrafteebene  entwickelt. 


1.  Introductory  Remarks 

The  fundamental  role  of  ground-combat  troops  (in  the  U.S.  Army's  own  words, 
e.g.  see  [l,  p.  i v ] )  is  to  "shoot,  move,  and  communicate."  Consequently  models 
of  combat  operations  must  in  some  manner  represent  the  attendant  processes  of 

3 

attrition,  movement,  and  C  (i.e.  command,  control,  and  communications) .  Tms 
paper  will  focus  on  the  modelling  of  combat  attrition,  although  some  consid¬ 
eration  does  have  to  be  given  to  the  other  two  processes  of  movement  and  C3, 
especially  as  they  influence  the  attrition  process.  More  specifically,  we  will 
examine  the  principal  methodologies  used  in  the  United  States  for  assessing 
casualties  in  simulated  combat  engagements,  with  emphasis  on  simulated  engage¬ 
ments  in  conventional  ground-combat  operations. 

The  two  attrition-modelling  approaches  that  we  will  examine  in  detail  are  as 
follows: 

(Al)  detailed  Lanchester- type^  models  of  attrition  in  tactical  engagements, 

(A2)  aggregated-force  casualty-assessment  models  based  on  the  use  of  index 
numbers  to  quantify  military  capabilities. 

Me  will  try  to  be  fairly  comprehensive  in  our  examination  of  these  attrition¬ 
modelling  approaches,  and  when  details  must  be  omitted,  references  to  further 
details  in  the  literature  will  be  given. 


^  So-called  after  the  pioneering  work  of  F.W.  Lanchester  [2).  We  will  refer 
to  any  differential-equation  model  of  combat  as  a  Lanchester-type  combat 
model  or  as  a  system  of  Lanchester-type  differential  equations  (or  some¬ 
times  simply  as  Lanchester-type  equations).  The  state  variables  are  typi¬ 
cally  the  numbers  of  the  various  different  weapon-system  types. 
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Although  scientific  work  on  the  analysis,  modelling,  and  evaluation  of  mili¬ 
tary  systems  in  simulated  combat  (see  Taylor  [3]  for  further  details)  has 
been  going  on  in  the  United  States  for  well  over  thirty  years  now,  there  un¬ 
fortunately  is  no  single  article  (or  even  small  number  of  introductory  arti¬ 
cles)  that  is  a  satisfactory  introduction  to  the  topic  of  attrition  modelling, 
especially  as  concerns  recent  developments  (and  there  have  been  some  very  im¬ 
portant  ones)  and  trends.  The  author  has  welcomed  the  opportunity  to  try  to 
integrate  and  synthesize  concepts  and  ideas  on  attrition  modelling.  Previous¬ 
ly,  such  ideas  and  results  have  been  widely  scattered  in  the  literature  (see, 
however,  survey  articles  by  Dolansky  (4)  and  Taylor  {3];  (5)). 

Combat  (especially  combat  between  company-sized  units  and  larger)  is  a  fan¬ 
tastically  complex  random  process.  This  inherent  complexity  of  the  combat  pro¬ 
cess  leads  to  great  complexity  in  operational  models  of  combat  attrition.  How¬ 
ever,  for  purposes  of  understanding  the  modelling  approaches  and  concepts  that 
may  be  used  to  build  such  operational  models,  it  is  convenient  to  abstract 
much  simpler  auxiliary  models  and  to  study  them.^  This  idea  of  using  simple 
auxiliary  models  to  illustrate  modelling  points  for  developing  and  understand¬ 
ing  complex  operational  models  apparently  has  not  been  clearly  articulated 
within  the  combat-modelling  community  (at  least  in  the  United  States),  and 
therefore  before  examining  currently  used  attrition-modelling  methodologies 
we  will  briefly  try  to  establish  a  conceptual  framework  for  reviewing  them. 


2.  Different  Types  of  Combat  Models 

Models  are  representations  and  idealizations  of  reality.  They  have  been  fairly 
widely  used  in  the  United  States  as  decision  aids  in  defense  planning  (e.g. 
to  evaluate  "on  paper"  proposed  weapon  systems  during  advanced  planning).  Fig¬ 
ure  1  shows  us  different  types  of  models  that  have  been  used  to  represent  com¬ 
bat  operations.  However,  for  present  purposes  let  us  focus  on  the  three  right¬ 
most  types  of  combat  models  shown  in  Figure  1  (see  Bonder  [7  )  for  further  de¬ 
tails  and  discussion): 

(Tl)  war  games, 

(T 2)  computer  simulations, 

(T3)  analytical  models. 

^  The  reverse  process  of  starting  with  a  simple  model  and  then  elaborating 
upon  it  and  enriching  it  in  details  is,  of  course,  the  approach  usually 
used  by  model  developers  to  build  their  models.  See  W.T.  Morris  (6  )  for  a 
lucid  discussion  of  this  enrichment  process. 
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Fig.  1:  The  Spectrum  of  Types  of  Combat  Models 


We  will  now  briefly  discuss  each  of  these  different  types  of  combat  models. 

The  distinguishing  feature  of  war  games  is  that  they  use  people  playing  roles 
to  simulate  decision  processes,  whereas  simulations  and  analytical  models  use 
algorithms  or  some  other  type  of  logic  to  represent  decision  processes.  Simu¬ 
lations  "act  out"  combat  operations,  usually  in  great  detail.  In  a  sense  they 
recreate  the  sand  table  (usually  with  the  help  of  a  digital  computer),  and 
battles  are  acted  out  on  this  automated  sand  table.  Such  combat  simulations 
usually  use  so-called  pseudo-random  numbers  to  determine  the  outcomes  of  ran¬ 
dom  events  (e.g.  the  outcome  of  firing  at  a  target)  and  are  consequently 
called  Monte  Carlo  simulations.  Analytical  models  use  symbols  to  represent 
the  combat  process  (and  hence  they  are  really  mathematical  models).  The  pro¬ 
cess  under  study  is  analyzed  and  abstracted  (i.e.  decomposed  into  basic  events 
and  activities).  Then  mathematical  submodels  of  events  and  activities  are  de¬ 
veloped  and  integrated  into  an  overall  structure.  Analytical  models  of  any 
degree  of  complexity  usually  do  not  yield  convenient  analytical  solutions  but 
require  numerical  approximation  methods  and  a  digital  computer  tor  the  gener¬ 
ation  of  numerical  results. 

Now  that  we  have  discussed  the  general  types  of  combat  models  that  are  avail¬ 
able  for  defense  analyses,  it  is  time  to  focus  on  the  types  of  models  that 
have  been  used  to  assess  outcomes  (in  particular,  casualties)  of  simulated 
tactical  engagements.  Thus,  the  three  approaches  that  are  currently  used  in 
the  United  States  for  assessing  casualties  in  simulation  tactical  engagements 
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are  as  follows  (see  Bonder  and  Farrell  [8]  for  further  details): 

(Al)  firepower  scores  (see  Stockfisch  [9,  pp.  6-27)), 

(A2)  Monte  Carlo  simulation  10  ;  11  , 

(A3)  analytical  models  (i.e.  differential  equations)  [8). 

The  firepower-score^  approach  is  an  index-number  method  for  aggregating  the 

heterogeneous  forces  on  a  side  into  a  single  equivalent  homogeneous  force.  We 

will  discuss  such  models  further  in  Section  6.1  below.  Lanchester-type  (i.e. 
4) 

differential-equation  ')  models  are  the  basic  type  of  analytical  model  used 
for  assessing  casualties,  although  other  types  of  analytical  submodels  may  be 
used  to  represent  processes  such  as  target  acquisition. 

If  we  try  to  relate  each  of  the  above  casualty-assessment  approaches  to  a  par¬ 
ticular  scale  of  combat  operations,  we  find  that  the  f i repower-score  approach 
and  Monte  Carlo  simulation  are  at  opposite  ends  of  the  spectrum  of  the  scale 
of  combat  operations  (i.e.  the  size  of  the  units  involved)  (see  Table  I).  The 
contents^  of  Table  I  are  only  generally  true,  with  exceptions  certainly  ex¬ 
isting.  As  we  see  from  this  table,  the  firepower  score  approach  has  been  pri¬ 
marily  used  for  casualty  assessment  in  large-scale  (i.e.  corps-level  and 


^  Indices  of  the  relative  combat  capabilities  of  military  units  (based  on  a 
"scoring  system"  for  the  weapons  employed  in  the  units)  have  been  used  by 
military  gamers  and  force  planners  in  the  United  States  for  at  least  thirty 
five  years.  We  are  here  generically  referring  to  both  such  indices  and  the 
associated  scores  as  firepower  scores.  (See  Stockfisch  [9)  and  below  in 
the  main  text  for  a  discussion  of  the  difference  in  meaning  between  the 
words  score  and  index  as  generally  used  in  defense  analyses.)  Members  of 
this  family  of  scores  and  indices  are  firepower  score/ index  of  combat  ef¬ 
fectiveness  (FS/ICE),  firepower  potential/unit  firepower  potential  (FP/ 
UFP),  firepower  potential  score/index  of  firepower  potential  (FPS/IFP), 
weapon  effectiveness  index/weighted  unit  value  (WFI/WUV),  weapon  effective¬ 
ness  value  (WEV),  antipotential  potential,  etc.  (see  Stockfisch  (9)  for 
further  references  and  a  guide  to  the  literature  about  firepower  scores; 
also  see  Honig  et  al.  ( 1 2 ,  Appendix  C  to  Chapter  III).  When  two  names  sep¬ 
arated  by  a  "slash")  are  given  above,  the  first  name  (e.g.  FS)  denotes  the 
scoring  system  for  weapon-system  types,  while  the  second  (e.g.  ICE),  iden¬ 
tifies  the  index  number  for  a  unit's  capability. 

4) 

'  We  are  calling  both  differential-equation  and  also  difference-equation  mod¬ 
els  Lanchester-type  models.  In  practice,  all  operational  models  of  combat 
systems  of  any  degree  of  complexity  use  finite-difference  methods  for  com¬ 
putation  and  thus  are  really  difference-equation  models.  However,  for  pur¬ 
poses  of  model  building,  it  is  much  more  convenient  to  think  in  terms  of 
differential  equations. 


(when  it  does  exist)  is  generally  poor.  The  following  documentation  and  in¬ 
formation  is,  however,  exceptionally  good  for  this  field.  General  informa¬ 
tion  about  contemporary  combat  models  in  the  United  States  is  available  in 
(14)  and  [ 1 5 ) .  Further  information  about  ATLAS  may  be  found  in  Kerlin  and 
Cole  (16)  or  (10),  while  that  about  CFf*  may  be  found  in  [17]  or  [  1 R ] .  Fur¬ 
ther  information  about  CAPf’ONETTI  may  be  found  in  / litmierman  [11]  nr  Adams 
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Table  I:  Combat-Assessment  Approach  Related 
to  Scale  of  Combat  Operations 


Modelling 

Approach 


Scale 

of 

Combat 


Example 


5) 


Model 


of 


firepower  score 

theater  -  ATLAS,  CEM 

Monte  Carlo 
Simultation 

infantry:  platoon  -  ASARS  II 

armor:  company/battalion- 
DYNTACS ,  CARMONETTE 

Lanchester-type 

Model 

battalion  -  BONDER/IUA 
division  -  DIVOPS 
theater  -  VECTOR- 2,  TWSP, 
DALFRAM,  DMEW 

See  footnote  on  previous  page. 


theater-level)  combat  models.  Although  there  are  exceptions,  high-resolution 
Monte  Carlo  simulation  has  been  a  feasible  assessment  approach  only  when  there 
have  been  no  more  than  about  100  elements  (e.g.  individual  tanks,  crew-served 
weapons,  etc.)  on  each  side.  On  the  other  hand,  Lanchester-type  models  have 
been  developed  for  the  entire  spectrum  of  combat  operations,  from  combat  bet¬ 
ween  company/battal ion-sized  units  to  theater-level  combat  operations. 

Thus,  we  find  that  Monte  Carlo  simulations  have  been  used  to  assess  casualties 
in  small-unit  combat  (i.e.  combat  between  battalion-sized  units  and  smaller), 
while  the  firepower-score  approach  applies  primarily  to  large-scale  (i.e. 
corps-level  and  theater-level)  combat.  However,  Lanchester-type  models  have 
been  developed  in  the  United  States  for  the  full  spectrum  of  combat  operations, 
from  small-unit  combat  to  large-scale  operations.  Thus,  if  one  wants  to  assess 
casualties  for  simulated  tactical  engagements  between  battalion-sized  units 
or  larger,  there  are  only  two  types  of  models  for  assessing  casualties  in  such 
tactical  engagements: 

(1)  detailed  Lanchester-type  models,  and 


5)  cntd. 

et  al.  (19),  while  that  about  DYNTACS  may  be  found  in  [ 20]  or  (21).  Infor¬ 
mation  about  BONDER/IUA  and  its  various  derivative  models  may  be  found  in 
(8);  (14);  (22);  (23);  (24),  while  that  about  DIVOPS  may  be  found  in  (25). 
The  theater-level  combat  model  named  VECTOR  is  documented  in  (26)  and  (27) 
DMEW  (see  (28))  is  also  a  theater-level  model,  as  is  TWSP  (see  Fain  et  al. 

(29))- 
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(2)  aggregated-force  models  based  on  quantifying  military  capabilities  with 
index  numbers  (i.e.  firepower-score  models). 


3.  Detailed  Models  versus  Aggregated-Force  Models  of  Attrition  in  Tactical 
Engagements 

Without  the  modern  high-speed  digital  computer  both  high-resolution  Monte 
Carlo  simulations  such  fis  OYNTACS  and  CARMONETTE  and  also  differential  combat 
models  such  as  BONDER/IUA  and  its  many  derivatives  would  be  impossible.  The 
modern  computer  provides  not  only  large-scale  memory  capacity  but  also  the 
ability  to  perform  millions  of  arithmetic  operations  per  second.  The  numeri¬ 
cal  integration  of  a  system  of  hundreds  of  ordinary  differential  equations  is 
consequently  possible.  Today  Lanchester-type  complex  system  models,  which  re¬ 
ly  on  modern  digital  computer  technology  for  their  implementation  (e.g.  see 
Bonder  and  Honig  [22]  or  [26]),  have  been  developed  for  various  levels  of  com¬ 
bat,  from  combat  between  battalion-sized  (e.g.  see  Bostwick  et  al.  [30]  or 
Hawkins  [23])  and  division-sized  [25]  units  to  theater-level  operations  (e.g. 
see  Cordesman  [31],  Farrell  [32],  or  [26];  [27]). 

All  the  above  complex  operational  models  that  are  conceptually  based  on  Lan¬ 
chester-type  equations,  however,  model  combat  attrition  in  detail  and  explic¬ 
itly  consider  the  many  different  weapon-system  types  that  can  be  individually 
attrited.  These  weapon-system  types  include  different  types  of  weapon  systems 
in  maneuver  units  and  different  types  of  fixed-wing  aircraft,  as  well  as  sep¬ 
arately  represented  field  artillery,  air  defense  artillery,  and  helicopter 
weapon  systems.  Such  Lanchester-type  models  represent  attrition  in  a  way  that 
reflects  the  internal  dynamics  of  combat  activities  and  relates  these  dynam¬ 
ics  to  specific  weapon-system  parameters  and  tactics  considered  important  in 
small-unit  engagements.  The  effects  of  individual  weapon-system  types  on  the 
outcome  of  a  theater-level  campaign  are  clearly  observable  and  bear  a  clear 
relationship  to  the  input  performance  assumed  (see  [27]  for  further  details). 
We  will  further  consider  such  models  in  Section  5  below. 

A  different  approach  for  modelling  attrition  in  large-scale  (i.e.  theater-lev¬ 
el  )  combat  operations  is  to  represent  attrition  in  a  macroscopic  fashion.  The 
many  different  weapon  systems  on  one  side  are  all  combined  together  into  a 
Single  scalar  quantity,  the  "combat  capability"  of  the  force,  and  combat  causes 
Attrition  of  this  index  number.  The  attrition  of  combat  capability  is  deter¬ 
mined  with  the  help  of  casualty-rate  curves  that  relate  the  relative  combat 
capabilities  of  the  forces  and  other  tactical  factors  to  their  casualty  rates 
(expressed  in  an  aggregated  fashion),  tosses  of  individual  weapon-system  types 
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are  then  determined  by  some  means  of  disaggregation.  Such  aggregated  loss- 
rate  relations  are  apparently  largely  judgmental ly  determined  (although  hav¬ 
ing  some  alleged  basis  in  empirical  combat  data),  and  the  author  knows  of  tio 
conceptual  approach  or  mathematical  models  for  relating  weapon-system-perfor¬ 
mance  parameters  and  other  operational  variables  to  the  numerical  determina¬ 
tion  of  these  aggregated-force  loss  rates. 


4.  Use  of  Simple  Auxiliary  Models  for  developing  and  understanding 
Complex  Operational  Models 

As  pointed  out  above  and  well  known  to  this  audience,  combat  (especially 
large-scale  combat  such  as  considered  for  NATO  studies,  e.g.  see  Huber  et  al. 
[33])  is  a  very  complex  process,  and  consequently  the  operational  models  that 
have  been  developed  in,  for  example,  the  United  States  are  also  very  complex. 
However,  for  purposes  of  discussing  and  understanding  such  complex  models  it 
is  much  more  convenient  to  consider  simplified  versions  of  them  than  to  con¬ 
sider  the  complex  models  themselves. 

Thus,  one  approach  for  understanding  the  reasons  why  a  large-scale  complex 
operational  model  produces  certain  output  results  for  particular  numerical 
input  data  is  to  abstract  a  simpler  model  (e.g.  one  with  fewer  variables  or 
simpler  functional  relations  between  them)  from  the  complex  one.  This  simple 
auxiliary  model  is  then  used  to  investigate  the  system  dynamics  of  the  more 
complex  model  by  considering  alternative  assumptions  and  data  estimates.  The 
simplified  auxiliary  model  should  be  intuitively  plausible  and  transparent 
but  yet  it  should  capture  the  basic  essence  of  the  complex  operational  model. 
This  idea  of  using  relatively  simple  auxiliary  models  in  conjunction  with  a 
complex  operational  model  is,  of  course,  not  new^ ,  but  the  author  knows  of 
no  clear  articulation  of  this  approach  for  understanding  large-scale  combat 
models. 


Geoffrion  [34]  has  suggested  a  similar  conceptual  approach  of  using  a  sim¬ 
ple  auxiliary  model  to  generate  tentative  hypotheses  to  be  tested  in  a 
full-scale  operational  model  and  thus  to  provide  guidance  for  further  (com¬ 
puterized)  higher-resolution  investigations.  We  also  have  felt  (Taylor 
[35])  that  the  use  of  relatively  simple  auxiliary  models  in  conjunction 
with  complex  operational  models  has  much  to  offer  for  the  analysis  of  mil¬ 
itary  operations  (see  also  Nolan  and  Sovereign  [36]  and  Weiss  [37]). 
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5.  Overview  of  Deteiled  Lenehester-Type  Models  of  Attrition  in  Tectical 
Engagements 


In  this  section  we  will  review  various  specific  models  and  associated  model¬ 
ling  methodology  for  detailed  Lanchester-type  models  of  attrition  in  tactical 
engagements.  Such  models  relate  casualties  sustained  to  both  the  tactical  and 
operational  variables  of  the  simulated  engagement  and  also  weapon-system  ca¬ 
pability  parameters.  We  will  illustrate  most  modelling  points  and  concepts  by 
considering  models  for  combat  between  two  homogeneous  forces.  Complex  opera¬ 
tional  models  that  have  been  developed  in  the  United  States  (see  Section  5.6 
below)  take  such  simple  models  as  their  point  of  departure  and  have  been  built 
by  the  process  of  enrichment  (see  Footnote  2)  and  Morris  [6]  for  further  de¬ 
tails). 

The  Lanchester-type  models  that  we  consider  here  are  all  deterministic  in  the 
sense  that  each  of  them  will  always  yield  the  same  output  for  a  given  set  of 
input  data.  Even  though  combat  between  two  military  forces  is  a  complex  ran¬ 
dom  process,  such  deterministic  combat  models  are  commonly  used  for  computa¬ 
tional  reasons  in  defense-planning  studies,  for  example,  to  assess  the  rela¬ 
tive  importance  of  various  weapon-system  and  force-level  parameters,  since 
they  give  essentially  the  same  results  for  the  mean  course  of  combat  as  do 
corresponding  stochastic  attrition  models. 


S.l  Various  Functional  Forms  for  Attrition  Rates 

Different  military/operational  situations  have  been  hypothesized  to  yield  dif¬ 
ferent  functional  forms  for  Lanchester-type  equations.  We  begin  by  reviewing 
two  simple  differential -equation  models  that  were  originally  considered  by 
F.W.  Lanchester  [2].  These  models  are  fundamental  for  representing  attrition 
and  still  may  be  considered  to  be  a  point  of  departure  for  modelling  combat 
attrition. 

Lanchester  [2]  hypothesized  that  under  conditions  "of  modern  warfare"  attri¬ 
tion  in  combat  between  two  homogeneous  forces  could  be  modelled  by 

g^  =  -ay  with  x (0 )  =  xQ  , 

S  (5.1.1) 

gf  =  -bx  with  y ( 0)  =  yQ  , 


with  x(t)  and  y(t)  denote  the  numbers  of  X  and  Y  at  time  t  after  the  battle 
begins,  and  a  and  b  are  positive  constants  that  are  today  called  Lanchester 
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attrition-rate  coefficients.  They  represent  the  effectiveness  of  each  side's 
fire,  i.e.  its  firepower.  We  will  call  (5.1.1)  Lanchester's  equations  for 
modern  warfare.  We  should,  of  course,  only  use  these  equations  when  both  x 
and  y  >  0  and,  for  example,  set  dx/dt  =  0  when  x  =  0.  Today  military  opera¬ 
tions  analysts  use  (5.1.1)  to  model  combat  attrition  under  the  following  two 
different  sets  of  circumstances: 

either  (a)  both  sides  use  "aimed"  fire  and  target-acquisition  times  are  con¬ 
stant,  independent  of  the  number  of  enemy  targets  (Weiss  (38]), 

or  (b)  both  sides  use  "area"  fire  and  a  constant-density  defense 
(Brackney  (39]). 

From  (5.1.1)  one  can  readily  deduce  Lanchester's  famous  "square  law" 

b(Xg  -  x2)  =  a(yg  -  y2)  ,  (6-1-2) 

which  yields  many  important  results,  e.g.  X  will  win  a  fight-to-the-f inish  if 
and  only  if  x^/y^  >  /a/b  (e.g.  see  Taylor  and  Comstock  (40 ] ) . 

Lanchester  also  hypothesized  that  when  both  sides  use  "area"  fire,  combat  at¬ 
trition  could  be  modelled  by 

^  =  -axy  ,  and  ^  =  -bxy  .  (5-1-3) 

It  should  be  noted  that  the  attrition-rate  coefficients  a  and  b  in  (5.1.1) 
and  (5.1.3)  represent  different  physical  quantities  and  are  consequently  giv¬ 
en  by  different  expressions  (see  Taylor  (41]).  For  want  of  a  better  alterna¬ 
tive,  we  will  call  (5.1.3)  Lanchester's  equations  for  area  fire,  even  though 
they  have  been  hypothesized  to  apply  to  certain  cases  of  "aimed"  fire  (see 
below).  Today  analysts  use  (5.1.3)  to  model  combat  attrition  under  the  fol¬ 
lowing  two  different  sets  of  circumstances: 

either  (a)  both  sides  use  "area"  fire  and  a  constant-area  defense  (Weiss 
(38],  Brackney  (39]), 

or  (b)  both  sides  use  "aimed"  fire  with  the  rate  of  target  acquisition 
being  inversely  proportional  to  the  number  of  enemy  targets  and 
also  being  the  controlling  (and  constraining)  factor  in  the  at¬ 
trition  process  (Brackney  (39]). 

From  (5.1.3)  one  can  also  deduce  Lanchester's  famous  "linear  law" 

b(*0  *  *)  =  a(yQ  -  y)  .  (5.1.4) 


which  also  yields  many  important  results. 
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Thus,  several  different  sets  of  physical  assumptions  may  be  hypothesized  to 
yield  the  same  functional  form  for  an  attrition  rate.  Consequently,  it  is 
more  convenient  to  refer  to  a  model  for  combat  between  two  homogeneous  forces 
in  terms  of  the  functional  forms  for  the  two  attrition  rates  than  to  refer  n 
terms  of  the  assumptions  as  we  have  done  above.  Let  us  therefore  introduce  a 
very  convenient  shorthand  for  referring  to  such  homogeneous-force  Lanchester- 
type  combat  models.  It  basically  involves  using  a  two-part  descriptor  X!Y, 
where  X  describes  the  attrition  rate  for  the  X  force  and  similarly  for  Y.  X 
and  Y  take  on  their  values  according  to  the  type  of  proportionality  for  the 
various  terms  in  a  side's  attrition  rate.  This  proportionality  is  expressed 
in  terms  of  the  number  of  firers  (denoted  as  F)  and/or  the  number  of  targets 
(denoted  as  T).  If  the  attrition  rate  is  independent  of  the  numbers  of  firers 
and  targets,  we  use  the  letter  C  (for  constant  attrition  rate).  When  there  is 
more  than  one  term  in  a  side's  attrition  rate,  the  same  approach  is  applied 
to  each  term,  with  a  plus  sign  separating  each  component  term  of  the  attri¬ 
tion  rate. 

Let  us  now  consider  some  examples  to  illustrate  this  shorthand.  For  example, 
for  Lanchester's  equations  of  modern  warfare  (5.1.1),  the  X  force's  attrition 
rate  is  (-dx/dt)  =  ay  so  that  it  is  proportional  to  only  the  number  of  enemy 
firers  (and  similarly  for  the  Y  force's  attrition  rate).  Consequently,  we  will 
refer  to  it  as  a  F|F  Lanchester-type  attrition  process  (or,  simply,  F ! F  at¬ 
trition).  Similarly,  Lanchester's  equations  for  area  fire  (5.1.3)  represent 
FT | FT  attrition,  since  each  side's  attrition  rate  is  proportional  to  the  pro¬ 
duct  of  the  numbers  of  firers  and  targets. 

Figure  2  shows  various  different  attrition-rate  functional  forms  that  have 
been  considered  in  the  literature  of  the  Lanchester  theory  of  combat.  We  have 
used  the  above  shorthand  notation  for  referring  to  these  various  attrition 
processes  in  the  figure.  Also  shown  for  each  process  are  the  state  equation 
(if  not  too  complicated)  and  the  first  person  (known  to  this  author)  to  have 
considered  it. 

Let  us  now  briefly  examine  various  sets  of  physical  assumptions  that  have 
been  hypothesized  to  yield  the  five  basic  attri tion-rate  functional  foff.t 
shown  in  Figure  2.  Conditions  hypothesized  to  yield  the  FjF  and  FT  1  FT  attr. 
tion  processes  have  been  discussed  above,  while  conditions  for  the  F!FT  pro¬ 
cess  (equivalently,  the  FT ) F  process)  are  just  a  combination  of  these  two 
sets,  with  one  set  applying  for  each  side.  For  example,  Brackney  39  has  hy¬ 
pothesized  that  the  F | FT  attrition  process  occurs  for  an  assault  by  the  X 
force  against  the  Y  force's  defensive  positions,  in  which  the  defenders  use 
aimed  fire  (with  X  targets  being  readily  acquired  by  virtue  of  their  "assault 
posture")  and  so  do  the  attackers,  only  their  sean  h  time  for  Y  tai  s  is 
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Fig.  2:  Various  functional  forms  for  attrition  rates 
that  have  been  considered  in  the  lanchester- 
combat- theory  literature 


relatively  large  (and  inversely  proportional  to  enemy  troop  density)  by  virtue 
of  the  enemy  remaining  under  cover  in  their  defensive  positions.  Also  Deitch- 
man  [42]  has  used  F | FT  attrition  to  model  insurgency  operations  (i.e.  guer¬ 
rilla  warfare)  in  which  Y-force  guerrillas  ambush  X-force  counterinsurgents. 
He  hypothesized  "aimed"  fire  for  the  ambushing  Y  force,  which  fires  on  the  X 
force  ("caught  in  the  open"),  but  that  the  ambushed  X  force  can  only  return 
area  fire,  since  its  members  do  not  know  the  exact  positions  of  individual  Y 
ambushers  and  consequently  return  fire  into  only  the  general  area  known  to  be 
occupied  by  the  enemy. 

Peterson  (43j;  (44]  has  hypothesized  that  T ( T  attrition,  i.e. 
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and  ^  •  (5.1.5) 

characterizes  the  early  stages  of  a  small-unit  engagement  in  which  the  vul¬ 
nerability  of  a  force  dominates  its  ability  to  acquire  enemy  targets.  He  in¬ 
troduced  this  model  to  extend  the  available  choice  of  basic  combat  models  and 
also  because  it  does  fit  limited  data  for  a  certain  type  of  engagement,  i.e. 
a  tactical  situation  in  which  all  weapons  of  the  two  forces  are  within  effec¬ 
tive  range  of  the  enemy  but  when  (due  to  cover,  concealment,  or  expert  camou¬ 
flage)  no  two  weapons  are  actually  intervisible.  Once  the  battle  actually  be¬ 
gins,  however,  this  model  no  longer  applies.  H.K.  Weiss  (45]  has  suggested 
that  force  vulnerability  may  become  the  dominant  factor  in  causing  losses  as 
combat  units  increase  in  size  and  become  increasingly  inefficient.  G.  Clark 
(46]  has  used  T|T  attrition  (5.1.5)  for  the  early  stages  of  a  small-unit  en¬ 
gagement  in  his  COMAN  model. 

The  last  attrition-rate  functional  form  shown  in  Figure  2  is  that  of  (F  +  T)| 
(F  +  T)  attrition,  i .e. 

=  -ay  -  fix  ,  and  ^  =  -bx  -  ay  .  (5.1.6) 

Two  situations  that  have  been  hypothesized  to  yield  the  above  equations  are 
(see  Figure  3): 

(51)  F | F  attrition  in  combat  between  two  homogeneous  forces  with  "operation¬ 
al"  losses  (Bach  et  al.  (47]  and  Morse  and  Kimball  (48]), 

(52)  F | F  attrition  in  combat  between  two  homogeneous  primary  forces  (see 
Weiss  (37])  with  superimposed  effects  of  supporting  fires  not  subject 
to  attrition  (Taylor  and  Parry  (49]). 

In  the  first  situation  (51),  for  example,  the  term  (fix)  in  X's  loss  rate,  i.e. 
(-dx/dt),  represents  "operational"  losses,  i.e.  losses  due  to  causes  other 
than  enemy  action  (47]  (e.g.  losses  due  to  sickness,  accidents,  desertions, 
etc.).  In  other  words,  the  model  holds  that  a  force  suffers  a  certain  amount 
of  casualties  because  of  its  very  size.  In  the  second  situation  (S2),  it  is 
assumed  that  F|F  attrition  occurs  between  the  primary  fighting  forces,  e.g. 
infantries,  and  that  the  supporting  weapons  employ  area  fire  against  enemy 
infantry  (again  see  Figure  3). 

A  general  form  for  homogeneous-force  attrition  rates  (and  which  yields  the 
square,  linear,  and  logarithmic  laws  as  special  cases)  has  been  given  by  P. 
Helmbold  (50],  who  proposed  a  modification  of  Lanrhester's  equations  for  "mod¬ 
ern  warfare"  to  account  for  inefficiencies  of  scale  for  the  larger  force  when 
force  sizes  are  grossly  unequal.  Helmbold  has  emphasized  that  Lanchester’s 
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Fig.  3:  Two  different  combat  situations  that  have  been 
hypothesized  to  yield  (F  ♦  T)|(F  +  T)  attrition 


classic  equations  for  modern  warfare  (5.1.1),  i.e.  the  F|F  attrition  model, 
imply  that  no  matter  how  unequal  the  force  levels  of  the  opposing  combatants 
may  be,  the  full  destructive  capability  of  each  side  can  always  be  focused  on 
the  enemy.  However,  sheer  limitations  on  available  spate  (to  say  nothing  of 
terrain-masking  and  reaction-time  effects)  may  well  prevent  the  larger  force 
from  using  its  full  destructive  capability  against  the  smaller  force,  espe¬ 
cially  if  the  opposing  force  is  much  smaller. 

Thus,  R.  Helmbold  (50 j  was  led  to  hypothesize  that  the  larger  force  suffers 
inefficiencies  of  scale  in  producing  casualties  when  force  sizes  are  grossly 
unequal.  He  confirmed  this  hypothesis^  by  analyzing  historical  data  (initial 


^  The  empirical  basis  for  Helmbold's  proposed  modification  is  not  reported 
in  (50j,  and  consequently  we  have  discussed  this  point  here.  Apparently  his 
CORG  reports  (51];  (52] ;  (53]  were  not  released  for  general  distribution 
to  the  public  until  after  the  publication  of  (50).  For  this  reason  we  have 
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and  final  force  levels  for  each  side)  for  several  hundred  battles  (see  Helm- 
bold  (5l)  ;  [52] ;  (53)).  Using  the  model  (5.1.1),  Helmbold  computed  the  com¬ 
batants'  relative  fire  effectiveness  a/b  for  each  battle  and  found  that  a/b 
and  the  initial  force  ratio  Xg/yg  were  strongly  positively  correlated  (see 
p.  7  of  (51)  and  pp.  31-35  and  58-59  of  [52]).  Thus,  as  the  initial  force  ra¬ 
tio  of  X  to  Y  Increases,  the  relative  fire  effectiveness  of  an  individual  X 
combatant  to  that  of  an  Individual  Y  one  decreases.  In  other  words,  a  side 
fights  less  "efficiently"  as  the  initial  force  ratio  (friendly  to  enemy 
forces)  increases. 

8ased  on  such  consideration  of  historical  combat  data,  Helmbold  [50)  intro¬ 
duced  a  modification  that  alters  relative  force-attrition  (or  fire-effective¬ 
ness)  capability  by  a  factor  depending  on  the  force  ratio.  For  temporal  vari¬ 
ations  in  fire  effectiveness,  his  proposed  modification  of  (5.1.1)  would  read 


jaf  =  -»(tKv(J)-y  *ith  x(°)  =  x0  • 

^  ^  -b(t)-Ex(*)-x  with  y(0)  =  yQ  , 


(5.1.7) 


where  a(t)  and  b(t)  denote  time-dependent  attrition-rate  coefficients,  and 
Ex  and  Ey  denote  the  fire-effectiveness-modification  factors  that  model  the 
inefficiencies  of  scale.  Helmbold  argued  that  the  fire-effectiveness-modifi¬ 
cation  factors  should  satisfy  the  following  three  requirements: 

(Rl)  Ex(u)  *  Ey(u)  =  E(u)  (i.e.  same  inefficiencies  of  scale  for  each  side), 
(R2)  E(u)  is  an  increasing  function  of  its  argument. 


(R3)  E( 1 )  =  1. 

Helmbold  then  considered  the  special  case  in  which  E(u)  is  a  power  function, 
i.e.  E(u)  =  uc  with  c  >  0.  In  this  case,  (5.1.7)  becomes 

^  =  -a(t)(J)1"W-y  withx(0)=x0, 

<  (5.1.8) 

%  -  -b(t)-(^)1'M-x  with  y ( 0 )  =  yQ  , 

where  we  will  call  W  the  "Weiss  parameter"  (see  [3)).  It  follows  that  W  =  1-c. 
We  will  refer  to  (5.1.8)  as  the  equations  for  Helmbold-type  combat. 


7)  entd. 

expounded  upon  how  Helmbold’ s  empirical  investigations  [ 5 1  ] ;  [52);  [53), 
which  were  all  done  well  before  the  submission  of  the  manuscript  [50),  mo¬ 
tivated  his  subsequent  proposed  modification  in  [ SO ) . 
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In  the  case  of  constant  attrition-rate  coefficients,  (5,1.8)  becomes 

dx  .  ,  x .  1  -W  ...  .... 

*  *d-(y)  -y  with  x(0)  =  Xq  , 

I  dy  t  ,y,l-W 


[5-1.9) 


with  y(Q)  =  yQ 


where  a  and  b  denote  constant  attrition-rate  coefficients.  Consequently,  we 

2  y.  l 

find  that  dx/dy  *  ( a/b ) (y/x )  ,  whence  the  state  equation  is  given  by 

b(*o“  -  *2w)  •  i(yZQ  -  y2w)  for  u  t  0  .  (5.1.10) 


b  jn(xn/x)  *  a  tn(yn/y)  for  W  *  0  . 


(5.1.11) 


Thus,  for  the  case  of  constant  attrition-rate  coefficients,  the  equations  for 
Helmbold-type  combat  yield  the  square  law  when  W  -  1,  the  linear  law  when 
W  *  1/2,  and  the  logarithmic  law  when  W  =  0.  Hence,  we  should  think  of  (5.1.9) 
as  a  general  combat  model  which  contains  many  of  the  classic  homogeneous- 
force  combat  models  as  special  cases. 


5.2  Fractional  Casualty  Rates  for  some  Selected  Models 

Many  aggregated-force  combat  models  assess  casualties  by  using  a  curve  (or 
table)  of  percent  daily  casualties  versus  the  force  ratio  of  the  attacker's 
combat  power  divided  by  that  of  the  defender  (with  other  tactical  variables 
held  constant).  In  Section  6.3  we  will  model  such  a  casual ty-rate  curve  with 
a  differential -equation  model.  Thus,  we  will  find  it  very  instructive  for  fu¬ 
ture  developments  in  this  paper  to  examine  casualty  rates  (expressed  as  a 
fraction  of  the  side's  current  strength)  for  several  of  the  above  simple  ho¬ 
mogeneous-force  attrition  models. 

Let  us  begin  by  examining  the  fractional  casualty  rate  for  a  side  in  combat 
modelled  by  Lanchester's  equations  for  modern  warfare  (5.1.1).  We  therefore 
consider,  for  example,  X ’ s  fractional  casualties  per  unit  time.  From  the  first 
of  equations  (5.1.1),  we  obtain 


(5.2.1) 


1  1  dx  \ 

/ X ' s  fractional  casualties  \ 

r  ;  a i)  1 

\  per  uni t  time  f 

where  u  denotes  the  force  ratio  of  X's  force  level  to  that  of  Y,  i.e.  u=  x/y, 
and  v  denotes  its  reciprocal,  i.e.  v  =  y/x. 

In  Figure  4  we  have  plotted  X's  fractional  casualties  per  unit  time  as  a 
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4: 


Relation  between  X's  casualty  rate  (expressed  as  a  fraction 
of  his  current  force  level  x(t))  and  the  force  ratio  (ex¬ 
pressed  as  the  ratio  of  the  attacker's  force  level  to  that 
of  the  defender)  for  Lanchester's  classic  model  of  "modern 
warfare"  dx 

3t  '  -a>- 


[NOTE:  In  the  legend  of  the  above  figure,  A  denotes  the 
attacker's  force  level,  and  D  denotes  that  of  the  defender.) 


function  of  a  certain  force  ratio.  The  force  ratio  that  we  have  used  is  the 
quotient  of  the  attacker's  strength  (here,  force  level)  divided  by  that  of 
the  defender  and  have  denoted  it  as  A/0,  since  most  combat  analyses  use  a 
force  ratio  of  this  type.  The  solid  line  in  Figure  4  represents  X's  fraction¬ 
al  casualties  per  unit  time  as  a  function  of  the  force  ratio  A/D  when  X  de¬ 
fends  and  Y  attacks.  It  is  a  straight  line  through  the  origin  with  a  slope 
equal  to  the  value  of  the  lanchester  attrition-rate  coefficient  a  as  the  read¬ 
er  can  see  by  referring  back  to  (5.2.1).  The  value  of  a  may  also  be  obtained 
from  this  figure  by  setting  v  -  1  in  (5.2.1)  and  reading  off  the  correspond¬ 
ing  value  for  X's  fractional  casualty  rate  from  the  curve.  Thus,  we  have  de¬ 
veloped  some  important  relations  between  the  fractional  casualty  rate  and  the 
Lanchester  attri tion-rate  coefficient.  Finally,  the  dashed  line  (which  is  a 
hyperbola)  in  Figure  4  represents  X's  fracMonal  rasualties  per  unit  time  as 
a  function  of  the  force  ratio  /*/ 0  in  the  other  case  in  which  X  attacks  and  Y 
defends.  Similar  curves  for  daily  casualty  rates  (hut  not  expressed  in  terms 
of  differential  equations)  are  commonly  used  to  assess  casualties  in  current¬ 
ly  operational  large-scale  aggregated- fori  e  combat  models  (see  Section  f). 
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Similarly,  for  the  Helmbold  type  equations  of  combat  (5.1.8)  we  have 


/  1  dx\ 

/X's  fractional  casualties) 

r  x  dt  /  ■ 

l  per  unit  time  J 

In  Figure  5  we  have  similarly  plotted  X's  fractional  casualties  per  unit  time 
versus  the  force  ratio  v  *  y/x  (denoted  in  the  figure  as  A/D)  for  the  case  in 
which  Y  attacks  and  X  defends.  In  this  figure  W  =  1  corresponds  to  the  case 
in  which  X's  casualty  rate  is  proportional  to  only  the  number  of  enemy  fir- 
ers,  and  (in  the  symmetric  case  in  which  Y‘s  casualty  rate  has  the  same  funct¬ 
ional  form)  consequently  the  corresponding  attrition  model  is  given  by  Lan- 
chester's  equations  for  modern  warfare  (5.1.1),  which  yield  the  square  law. 

As  before,  we  see  that  in  this  case  (i.e.  when  W  *  1)  X's  fractional  casual¬ 
ties  per  unit  time  are  directly  proportional  to  the  force  ratio  A/D  when  Y 
attacks  and  X  defends.  Referring  back  to  the  first  of  equations  (5.1.9),  we 
see  that  M  =  Wj  corresponds  to  a  more  efficient  use  of  firepower  for  force 
ratios  v  *  A/D  *  y/x  >  1  than  does  W  -  U2  when  1  2  Wj  >  W2>  since  the  corre¬ 
sponding  fire-effectiveness-modification  factor  for  W  =  W,  (i.e.  Ew(x/y)  = 

1-Wi  1  " 

(x/y)  A)  is  greater  than  that  for  W  =  W2  when  y/x  >  1.  Figure  6  shows  the 

same  type  of  plot  when  X  is  the  attacker  and  Y  is  the  defender.  In  this  case, 

the  casualty-rate  curve  corresponding  to  the  square  law  is  a  hyperbola. 


fORCE  R*no,  */0 


Fig.  5:  Relation  between  X’s  fractional  casualty 
rate  and  the  force  ratio  for  the  model 
dx  x  1 -W 

*  -a-(-)  y  when  X  defends 


FORCE  RATIO,  A/D 

Fig.  6:  Relation  between  X's  fractional  casualty 
rate  and  the  force  ratio  for  the  model 

^  r  when  X  attacks 

S.3  Determination  of  Numerical  Values  for  Attrition-Rate  Coefficients 

For  applying  any  kind  of  detailed  Lanchester-type  attrition  model  to  study  a 
particular  hypothesized  combat  engagement  in  a  defense-planning  study,  one 
must  be  able  to  predict  rates  at  which  weapon  systems  would  inflict  and  sus¬ 
tain  casualties.  Within  the  context  of  an  assumed  functional  form  for  the  at¬ 
trition  rates,  this  means  determining  numerical  values  for  the  associated 
Lanchester  attrition-rate  coefficients.  Two  approaches  that  have  been  used  in 
the  United  States  for  such  numerical  determinations  are  based  on  using 

(Al)  a  statistical  estimate  based  on  "combat"  data  generated  by  a  detailed 
Monte  Carlo  combat  simulation,  and 

(A2)  an  analytical  submodel  of  the  attrition  process  for  the  particular  com¬ 
bination  of  firer  and  tarqet  types. 

The  first  approach  (Al)  has  been  called  by  Bonder  (7]  the  use  of  a  fitted-pa¬ 
rameter  analytical  mode),  since  the  basic  idea  is  to  statistically  estimate 
parameters  for  the  attrition-rate  coefficient  from  the  output  of  a  high-reso¬ 
lution  Monte  Carlo  combat  simulation  (see  F igure  7).  The  "combat"  data  or 
outputs  of  the  simulation  are  used  to  fit  one  or  more  free  parameters  in  the 
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Fig.  7:  Basic  idea  of  the  fitted-parameter  analytical  model 
for  complimentary  use  of  Monte  Carlo  simulation  and 
Lanchester-type  attrition  models 


analytical  model  so  that  it  will  at  least  duplicate  and  hopefully  predict  re¬ 
sults  comparable  to  those  obtainable  from  the  simulation  model.  The  COMAN 
model  (see  Clark  [ 46 1 )  is  an  example  of  such  a  fitted -parameter  model.  En¬ 
couraging  results  have  been  reported  (see  [10]). 

S.  Bonder  [ 7  ]  has  called  the  second  approach  (A2)  the  use  of  a  freestanding 
or  independent  analytical  model,  since  this  type  of  analytical  model  can  be 
run  independently  of  any  detailed  Monte  Carlo  simulation  of  the  combat  pro¬ 
cess.  The  basic  conceptual  idea  is  to  develop  an  analytical  expression  for 
each  attrition-rate  coefficient  by  considering  a  single  firer  engaging  a "pas¬ 
sive"  target  (i.e.  one  that  does  not  fire  back).  One  designs  such  a  model  to 
use  the  same  types  of  inputs  as  used  by  detailed  Monte  Carlo  simulations  of 
the  same  combat  process.  An  example  of  such  an  independent  analytical  model 
is  the  BONDER/IUA  differential  combat  model,  which  was  first  used  in  the  Unit¬ 
ed  States  in  1969  [7],  and  the  many  subsequently  enriched  versions  of  it  (see 
Footnote  5)  for  references). 

Bonder  and  Farrell  [8]  have  used  the  second  approach  (A2)  to  develop  general 
methodology  for  determining  attrition-rate  coefficients  for  a  wide  spectrum 
of  weapon-system  types  engaging  specified  target  types.  Basically,  their  ap¬ 
proach  is  founded  upon  calculation  of  a  Lanchester  attrition-rate  coefficient 
as  being  the  reciprocal  of  the  expected  time  for  an  individual  firer  to  kill 
a  single  target  (see  Barfoot  [54]).  For  the  model  (5.1.1)  this  means 

a  *  i  ,  (5.3.1) 

M  xyl 

where  TXy,  a  random  variable  (frequently  abbreviated  as  r.v.),  denotes  the 
time  required  for  an  individual  Y  firer  to  kill  a  single  X  target,  and  E[T] 
denotes  the  expected  value  of  T  (a  r.v.). 
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To  facilitate  such  analysis  Bonder  and  Farrell  [B]  have  classified  the  en¬ 
gagement  of  particular  target  types  by  different  weapon-system  types  accord¬ 
ing  to  the  taxonomy  shown  In  Table  II.  Weapon-system  types  are  first  classi¬ 
fied  according  to  the  mechanism  by  which  they  kill  particular  target  types 
(l.e.  their  lethality  characteristics)  as  being  either  Impact-to-klll  systems 
or  area-lethality  systems.  Within  each  of  these  two  categories  Bonder  and 
Farrell  further  classify  weapon-system  types  according  to  how  they  use  firing 
Information  to  control  the  system's  aim  point  and  their  delivery  characteris¬ 
tics,  l.e.  the  firing  doctrine  employed.  Expressions  have  been  developed  for 
Lanchester  attrition-rate  coefficients  corresponding  to  the  weapon-system 
classifications  tagged  with  asterisks  *  in  Table  II. 


Table  II:  Classification  of  Weapon-System  Types  for 

the  Development  of  Lanchester  Attrition-Rate 
Coefficients  for  the  Model  (5.1.1) 

(from  Bonder  and  Farrell  [8]) 


Lethality  Mechanisr 

(1) 

Impact 

(2) 

Area 

Firing  Doctrine 

a) 

Repeated  Single  Shot 

(a)* 

Without  Feedback  Control  of  Ain  Point 

(b)* 

With  Feedback  on  Innediately  Preceding  Round 
(Markov-Dependent  Fire) 

(c) 

with  Complex  Feedback 

(2) 

Buret 

Fire 

• 

(a) 

Without  Ain  Change  or  Drift  in  or  Between  Bursts 

<b>* 

tfith  Aim  Drift  in  Burets,  Aim  Refixed  to  Original 
Aim  Point  for  Each  Burat 

(c) 

With  Ain  Drift,  Re-ain  Between  Bursts 

(3) 

Multiple  Tube  Firing:  Feedback  Situationa  (la), 

(lb),  (lc) 

(el* 

Salvo  or  Volley 

(41 

Mixed 

-Mode  Firing 

(a) 

Adjustment  Followed  by  Multiple  Tube  Fire 

(bl* 

Adjustment  Followed  by  Burst  Fire. 

'indicates  that 
by  Bonder  and 

analysis  of  this  category  has  been  performed 

Farrell  (01. 
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Let  us  now  illustrate  Bonder's  attrition-rate-coefficient  methodology  for  the 
model  (5.1.1).  A  large  class  of  weapon  systems  (e.g.  tanks  firing  at  tanks, 
anti-tank  weapon  systems  firing  at  tanks,  etc.)  may  be  classified  as  Markov- 
dependent-  fire  weapons,  i.e.  the  outcome  of  the  firing  of  a  round  by  the  weap¬ 
on  system  depends  on  only  the  outcome  of  the  immediately  preceding  round.  For 
such  weapon  systems  and  an  impact-to-kill  lethality  machanism,  Bonder  ( 55 ] ; 
[56]  (see  also  Kimbleton  [57])  has  developed  the  following  expression  for  the 
expected  time  for  an  individual  firer  to  kill  an  enemy  target 

ElTl  =  +  ei  '  *h  4  PfKlHy 

+  PHilmP  {”  '  +  p<hlh>  ‘  <5-3-2) 

where  all  variables  are  defined  in  Table  III.  The  corresponding  Lanchester 
attrition-rate  coefficient  is  then  given  by  the  reciprocal  of  (5.3.2),  e.g. 
see  (5.3.1)  above.  The  above  ideas  are  the  conceptual  basis  for  all  attrition 
calculations  in  the  BONDER/IUA  model.  Let  us  finally  record  here  that  the  pre¬ 
cise  conditions  for  (5.3.2)  to  hold  are  essentzilly  given  by  the  following 
assumptions: 

(Al)  Markov-dependent  fire  with  parameters  p,  P(hlh),  and  P(hlm), 

(A2)  geometric  distribution  for  the  number  of  hits  required  for  a  kill  with 
parameter  P(KIH). 

Although  (5.3.2)  is  a  rather  formidable-looking  expression,  it  is  easily  eval¬ 
uated  with  an  automated  algorithm  such  as  is  available  for  a  computerized  mod¬ 
el  such  as  the  BONDER/IUA.  Moreover,  in  several  special  cases  of  interest 
(5.3.2)  simplifies  appreciably.  In  the  simplest  case,  we  have  (1)  target-ac¬ 
quisition  time  negligible  (i.e.  set  ta  *  0),  (2)  uniform  rate  of  fire  (i.e. 
ti  *  t^  *  tm  =  l/v,  where  v  is  simply  the  firing  rate),  (3)  statistical  inde¬ 
pendence  among  firing  outcomes  (i.e.  p  =  P(hlh)  =  P(h|m)  =  P^SH),  and  (4) neg¬ 
ligible  time-of-flight  of  the  projectile  (i.e.  set  tf  =  0);  and  then  (5.3.2) 
reduces  to 

E(T]  =  l/(v  PSSK),  (5.3.3) 

where  the  single-shot  kill  probability  P^x  is  given  by  P^  =  pssh'p(k,h)- 
In  this  case,  the  Lanchester  attrition-rate  coefficient  is,  for  example,  con¬ 
sequently  given  by 

a  *  vypsskxy* 


(5.3.4) 
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Table  III:  Factors  Included  In  Expression  for  Lanchester 
Attrition-Rate  Coefficient  for  Single-Shot 
Markov-Dependent-Fire  Weapon  Systems  with  a 
Geometric  Distribution  for  the  Number  of  Hits 
Required  for  a  Kill 


Tlm«  to  acquire  a  target,  t# 

Time  to  fire  flrat  round  after  target  acquired,  t^ 

Time  to  fire  a  round  following  a  hit,  t ^ 

Time  10  fire  a  round  following  a  mlaa,  tffl 

Time  of  flight  of  the  projectile,  tf 

Probability  of  a  hit  on  flrat  round,  p 

Probability  of  a  hit  on  a  round  following  a  hit,  P(h|h) 

Probability  of  a  hit  on  a  round  following  a  mist,  P(h!m) 

Probability  of  destroying  a  target  given  it  la  hit,  P(k|H) 


which  is  a  very  intuitively  appealing  result. 

It  is  very  instructive  to  combine  an  idea  due  to  Brackney  [39]  concerning  the 
target-acquisition  process  with  Bonder's  general  expression  (5.3.2).  In  this 
case,  we  assume  that  (1)  the  mean  time  to  acquire  a  target  is  inversely  pro¬ 
portional  (let  k  denote  the  constant  of  proportionality)  to  target  density, 
and  also  (2)  through  (4)  above  in  the  previous  paragraph.  Then,  for  example, 
the  expected  time  to  kill  a  target  reduces  to 


ElTsvl 


kY*X  +  1 

X  vYPSSKxy  ’ 


(5.3.5) 


where  Ax  denotes  the  "presented"  area  occupied  by  the  X  force  and  which  is 
visually  searched  by  a  Y  firer.  When  the  first  term  predominates  (i.e.  the 
constraining  factor  in  killing  targets  is  acquiring  them),  we  may  set 
E[Txy]  *  kYAx/x  and  obtain 


dx  - 

it  =  “axy  •  (5.3.6) 

where  a  =  l/(kyAx).  Thus  (as  pointed  out  by  Brackney  (39)),  the  assault  of  a 
defensive  position  may  be  more  appropriately  modelled  with  the  equations  of  a 
FT | F  attrition  process. 

As  seen  In  (5.3.1)  and  (5.3.2),  the  Lanchester  attrition-rate  coefficient  de¬ 
pends  directly  on  the  factors  given  in  Table  III.  However,  it  also  depends 
Indirectly  on  variables  such  as  the  range  between  firer  and  target,  target 
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posture,  terrain,  target  motion,  etc.  Many  people  (e.g.  Bonder  and  Farrell 
(8])  feel  that  for  many  tactical  situations  the  principal  factor  is  the  range 
between  firer  and  target.  For  illustrative  purposes,  let  us  therefore  consid¬ 
er  that  the  Lanchester  attrition-rate  coefficients  in  (5.1.1)  explicitly  de¬ 
pend  on  range, i.e. 


a  =  a(r),  and  b  =  8(r), 


(5.3.7) 


where  r  denotes  the  range  between  firer  and  target.  Bonder  and  Farrell  [8, 
pp.  196-200]  have  used  expressions  like  (5.3.1)  and  (5.3.2)  to  examine  the 
variations  in  weapon-system  kill  rate  with  range.  The  following  functional 
form  fits  data  for  a  number  of  representatives  weapon  systems  (see  Bonder  and 
Farrell  18]  for  further  details) 


a(r)  = 


for  0  £  r  £  r  , 
a 

for  r  £  r  , 
a 


(5.3.8) 


where  ra  denotes  the  maximum  effective  range  of  Y's  weapon  system,  aQ  >  0, 
and  u  £  0.  The  constant  w  allows  us  to  model  the  range  dependence  of  the  weap¬ 
on  system's  kill  rate  (see  Figure  8). 

When  a  weapon  system  employs  "area"  fire  and  enemy  targets  defend  a  constant 
area,  the  expression  for  the  Lanchester  attrition-rate  coefficient  takes  a 
different  form  and  also  depends  (among  other  things)  on  the  vulnerable  area 
of  the  target  (denoted  as  a v)  and  the  lethal  area  of  the  projectile  fired  by 
the  firer' s  weapon  system  (denoted  as  aL).  In  general  a  rather  complicated 
expression  is  obtained  for  such  an  attrition-rate  coefficient  (e.g.  see  Bonder 
and  Farrell  [8,  pp.  141-162]),  but  in  special  cases  these  simplify  (cf.  equa¬ 
tion  (5.3.4)  above),  e.g.  for  "small  arms  fire"  when  av  »  aL  and  for  a "weap¬ 
on  of  great  lethality"  when  aL  »  ay. 

Thus,  two  cases  in  which  simple  expressions  are  obtained  for  attrition-rate 
coefficients  for  "area"  fire  and  a  constant-area  defense  are  as  follows  for: 


(1)  small  arms  fire  (i.e.  ay  »  aL),  and 

(2)  weapons  of  large  lethality  (i.e.  aL  »  av). 


)  In  actual  application,  a  model  like  BONDER/IUA  or  VECTOR-1  allows  such  at¬ 
trition-rate  coefficients  to  additionally  depend  on  other  variables  (e.g. 
target  posture,  terrain,  etc.)  that  depend  on  the  positions  of  firers  and 
targets  and  other  operational  factors  and  that  may  change  over  time. 
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Figure  8:  Dependence  of  Y‘s  attrition-rate  coefficient  a(r)  on  the 
exponent  y  with  the  maximum  effective  range  of  the  weapon 
system  and  kill  rate  at  zero  range  held  constant. 

(NOTES:  1.  The  maximum  effective  range  of  the  system  is 
denoted  as  r  =  2000  meters.  2.  a(0)  =  oq  =  0.6  X  casual- 
ties/(unit  time  v  number  of  Y  firers)  denotes  the  weapon- 
system  kill  rate  for  Y  at  zero  force  separation  (range). 
3.  The  opening  range  of  battle  (see  Section  6.4)  is  de¬ 
noted  as  rg  *  1250  meters  and  (as  shown)  rg  <  ra.  ) 


In  both  cases,  we  may  consider  the  assumptions  of  "area"  fire  and  a  constant 
area  defense  by  targets  to  yield 


dx  p 

M  =  "VY  *SSKxy 


(5.3.9) 


For  small  arms  fire,  we  may  calculate  pssKXy  by  consider  a  lethal  "dot"  being 
randomly  placed  into  a  large  region  (of  area  Ax)  that  contains  x  "vulnerable 
circles"  (each  of  area  ayx)  randomly  placed  in  the  region.  It  follows  that 

**SSK  «  ay  x/Ax  and  thus 


dx 


vyav 


(5.3.10) 


For  weapons  of  large  lethality  (e.g.  artillery),  we  may  calculate  pSSKxy  by 
considering  a  "lethal  circle"  being  placed  into  the  region  which  contains  x 
"vulnerable  dots".  If  there  are  no  multiple  kills,  then  *  a^  s/Ax  and 
thus 


XY 


dx 

11  m 


xy, 


(5.3.11) 


j 


164 


James  G.  Taylor 


where  denotes  the  lethal  area  of  the  projectile  fired  by  Y‘s  weapon  sys¬ 
tem.  A  more  refined  analysis  is  contained  in  Bonder  and  Farrell  [8]  (see  also 
Cherry  [58]). 


5.4  Additional  Operational  Factors  to  be  considered  in  Attrition  Models 

In  order  to  reflect  more  of  the  complexity  of  combat  we  can  enrich  the  above 
simple  attrition  models  by  considering  additional  operational  factors  such 
as: 

(FI)  range-dependent  weapon-system  capabilities, 

(F2)  other  temporal  variations  in  fire  effectiveness, 

(F3)  target-acquisition  considerations, 

(F4)  command,  control,  and  communications, 

(F5)  unit  breakpoints, 

(F6)  unit  deterioration  due  to  attrition, 

(F7)  suppressive  effects  of  weapon  systems, 

(F8)  effects  of  logistics  constraints. 


The  addition  of  such  factors  to  attrition  models  greatly  enhances  their  oper¬ 
ational  realism.  Although  we  will  not  in  all  cases  be  able  to  report  full  de¬ 
tails  of  computational  results,  it  turns  out  that  such  enriched  models  usual¬ 
ly  yield  quite  different  results  than,  for  example,  Lanchester’s  classic  con¬ 
stant-coefficient  equations  for  a  F|F  attrition  process. 

S.  Bonder  [59];  [60]  (see  also  Bonder  and  Farrell  [8])  has  stressed  the  im¬ 
portance  for  evaluating  many  types  of  proposed  weapon  systems  of  using  vari¬ 
able-coefficient  combat  models  to  represent  temporal  variations  in  firepower 
on  the  battlefield.  Such  a  case  occurs,  for  example,  when  the  range  between 
firers  and  targets  changes  appreciably  during  battle.  We  have  already  pointed 
out  above  (see  Section  5.3)  that  attrition-rate  coefficients  should  generally 
be  considered  to  depend  on  firer- target  range. 

Let  us  therefore  consider  "aimed-fire"  combat  between  two  homogeneous  forces 
and  assume  that  target-acquisition  times  do  not  depend  on  the  numbers  of  tar¬ 
gets.  We  further  assume  that  one  force  attacks  at  constant  speed  the  other 
force's  static  defensive  position.  Then  (for  x  and  y  >  0)  we  have 


dx 

m  = 


-a(r)y 


< 


with  x(t=0)  =  xQ  , 


-8(r)x 


with  y(t=0)  =  yQ  , 


(5.4.1) 
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where  r  denotes  the  range  between  opposing  forces  (see  Figure  9),  and  a(r) 
and  e(r)  denote  range-dependent  attrition-rate  coefficients.  Range  is  relat¬ 
ed  to  time  by 

r(t)  =  rQ  -  vt  ,  (5.4.2) 

where  rQ  denotes  the  opening  range  of  battle  and  v  >  0  denotes  the  constant 
attack  speed. 


Figure  9:  Diagram  of  Bonder's  constant-speed  attack  model. 

Force  separation,  r(t),  is  given  by  r(t)  =  r0  -  vt. 


Let  us  now  illustrate  the  importance  of  such  range-dependent  weapon-system 
capabilities.  We  assume  that  the  Y  weapon  system’s  kill  rate  a(r)  is  given  by 
(5.3.8)  and  similarly  for  8(r)  with  parameters  rg  and  v.  In  Figure  10  we  have 
plotted  force-level  trajectories  for  three  different  battles,  denoted  as  (A), 
(B)  and  (C).  In  all  these  battles  both  weapon  systems  have  the  same  maximum 
effective  range,  i.e.  r^  =  rg  =  rg,  and  the  battle  begins  at  this  range,  i.e. 

rn  =  r  .  For  these  battles  we  have  held  the  kill  rates  at  zero  force  separa- 

o  e 

tion,  i.e.  Og  =  a(0)  and  fig,  constant  and  have  varied  between  battles  the 

manner  in  which  a(r)  and  e(r)  depend  on  range,  i.e.  for  0  s  r  <  rf 

(A)  constant-constant:  a(r)  =  aQ  and  e(r)  =  Bg  , 

(B)  linear-linear:  a(r)  =  aQ(l-r/re)  and  a( r )  =  8g(l-r/re)  , 

(C)  linear-quadratic:  a (r)  =  ag(l-r/rp)  and  B(r)  =  Bg(l-r/re)2  . 

We  see  from  Figure  10  that  battle  outcome  may  be  quite  sensitive  to  the 
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Figure  10:  Force-level  trajectories  of  X  and  Y  forces  for  three 
different  battles  denoted  in  the  figure  as  (A),  (B) 
and  (C)  and  explained  in  the  main  text  with  each 
side's  fire  effectiveness  modelled  by  the  power  at¬ 
trition-rate  coefficients  for  rg  =  ra  *  re  »  re  = 

2000  meters,  ag  *  0.06  X  (casualties/minute)  per  Y 
firer,  gg  »  0.6  Y  (casualties/minute)  per  X  firer, 
v  *  5  mpn,  xg  *  10,  and  yg  *  30.  The  symbol  *  denotes 
the  end  of  a  force-level  trajectory  due  to  annihilation 
of  the  enemy  force. 


variation  in  weapon- system  kill  rate  with  range.  This  type  of  battle  was  first 
studied  by  Bonder  (59];  [60],  and  it  also  appears  as  an  example  in  much  of  the 
author's  work  (see  Taylor  [3];  [41];  161],  Taylor  and  Brown  [62];  [63],  and 
Taylor  and  Comstock  [40]). 

If  we  use  (5.4.2)  to  eliminate  range  from  the  attrition-rate  coefficients  in 
(5.4.1)  we  obtain  the  following  model  with  time  as  the  independent  variable 
(see  Taylor  [3];  [41]) 

• 

S  -  *a(t)y 
W  -  -b(t)x 


with  x(0)  =  x0  , 
with  y (0)  *  y0  , 


(5.4.3) 


where  the  time-dependent  attrition-rate  coefficients  are  given  by 
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a(t)  •  k#(t  ♦  C)v,  and  b(t)  *  kfa(t  +  C  +  D)u,  (5.4.4) 


with  C  ■  (ra  -  Tq)/v  and  D  *  (rg  -  ro)/v.  Moreover,  we  should  use  a  model 
like  (5.4.3)  instead  of  (5.1.1)  when  there  are  temporal  variations  in  fire¬ 
power  due  to  changes  over  time  in  firing  rates,  target-acquisition  rates,  tar¬ 
get  postures  and  hence  vulnerabilities,  etc.  [i.e.  the  inputs  for  calculating 
the  attrition-rate  coefficients  change  over  time;  see  (5.3.1)  and  (5.3.2)). 

An  operational  model  such  as  BONDER/ 1UA  actually  recomputes  such  attrition- 
rate  coefficients  at  every  "time  step"  (used  for  numerical  integration  of  the 
system  of  differential  equations  for  casualty  assessment)  during  the  battle. 

Another  important  factor  to  consider  is  engagement  termination,  i.e.  when  to 
break  off  the  engagement  and  stop  casualties.  Although  there  is  some  conflict¬ 
ing  evidence  (see  Helmbold  (64)),  most  military  planners  in  the  United  States 
usually  assume  that  a  ground-combat  unit  will  break  off  an  engagement  and  try 
to  disengage  from  the  enemy  when  it  has  suffered  a  certain  percentage  of  ca¬ 
sualties  for  a  company-sized  unit  in  the  attack  and  fifty  percent  in  the  de¬ 
fense).  Thus,  the  simplest  model  of  engagement  termination  is  to  assume  that 
battle  outcome  depends  deterministically  on  only  the  force  levels:  the  first 
side  to  be  reduced  to  a  previously  specified  force  level  (his  "breakpoint") 
is  assumed  to  break  off  the  engagement  (he  "loses").  Incorporating  this  idea 
into  the  simple  model  (5.1.1),  we  obtain 


for  x  >  Xgp  and  y  >  yBp, 
otherwise, 

for  x  >  xBp  and  y  >  yBp, 
otherwise. 


(5.4.5) 


where  xBp  denotes  X's  breakpoint  force  level  and  similarly  for  yBp. 

The  inclusion  of  unit  "breakpoints"  in  attrition  models  is  very  important, 
since  battle  outcome  (i.e.  who  "wins"  and  who  "loses")  is  very  sensitive  to 
the  choice  of  breakpoints.  To  Illustrate  this  point,  let  us  observe  that  Y 
will  “win"  an  engagement  modelled  by  (5.4.5)  if  and  only  if 


N) 

*0 


a 

F 


(5.4.6) 
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X  Y 

where  xBp  =  fgpx0  and  ^BP  *s  s’m^ar1>'  defined.  Table  IV  shows  us  the  sensi¬ 
tivity  of  battle  outcome  to  the  selection  of  breakpoints.  Clearly,  more  re¬ 
search  (especially  on  empirical  considerations)  is  necessary  on  various  as¬ 
pects  of  modelling  the  engagement-termination  process. 


Table  IV:  Influence  of  Unit  Breakpoints  on  the  Outcome  of  Battle  for 
an  Attack  by  X  Against  V  with  Battle  Dynamics  Given  by 
Lanchester's  Equations  for  Modern  Warfare 
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As  we  saw  in  Section  5.3,  target  acquisition  is  reflected  in  the  simple  model 
for  the  Lanchester  attrition-rate  coefficient  through  ta,  which  appears  in 
(5.3.1)  through  (5.3.2)  see  also  the  "Brackney  model"  (5.3.5)  .  However, 
Vector  Research,  Inc.  (see  (26,  pp.  103-108]  or  [27,  pp.  43-45])  has  devel¬ 
oped  a  more  refined  (i.e.  enriched  in  operational  details)  model  for  the  tar¬ 
get-acquisition  process  and  its  Impact  on  attrition-rate  coefficients.  They 
consider  that  the  two  major  factors  determining  the  value  of  an  attrition- 
rate  coefficient  are  (1)  the  acquisition  and  selection  of  targets,  and  (2) 
the  conditional  kill  rate  (i.e.  the  rate  at  which  acquired  targets  are  de¬ 
stroyed).  Concerning  target  acquisition  and  selection,  the  proportion  of  time 
that  a  weapon  is  actively  engaging  an  enemy  target  depends  on  the  interaction 
of  three  processes: 

(PI)  the  line-of-sight  process  (which  determines  when  a  given  target  is  vis¬ 
ible  or  invisible  to  a  potential  firer), 

(P2)  the  target-acquisition  process  (which  determines  the  time  required  for 
a  firer  to  acquire  a  particular  target),  and 

(P3)  the  target  selection  process  (which  specifies  a  scheme  by  which  a  weap¬ 
on  crew  chooses  to  engage  a  particular  target  from  among  those  that 
have  been  acquired). 

The  exact  way  in  which  the  above  three  processes  interact  depends  in  an 
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essential  way  on  which  of  two  kinds  of  acquisition  and  target-selection  modes 
the  weapon  systems  employ  -  serial  or  parallel  acquisition  (see  (26)  or  (27) 
for  further  details).  Suppressive  effects  of  weapon  systems  may  be  accommo¬ 
dated  in  Vector  Research's  model,  although  the  phenomenological  basis  of  sup¬ 
pressive  effects  is  poorly  understood  at  this  time  (see  [65]). 


S.S  Modelling  Attrition  for  Combat  between  Heterogeneous  Forces 

We  have  considered  above  various  aspects  of  attrition  modelling  for  combat 
between  two  homogeneous  forces,  but  actual  combat  consists  of  many  different 
weapon-system  types  operating  together  as  "combined-arms  teams".  For  example, 
there  may  be  infantry  (armed  with  several  types  of  weapons),  tanks,  artillery, 
mortars,  etc.  on  each  side.  Let  us  therefore  consider  combat  between  such  het¬ 
erogeneous  forces  and  briefly  indicate  how  the  above  basic  ideas  on  modelling 
combat  attrition  are  extended  and  adapted  to  such  cases. 

For  illustrative  purposes,  we  consider  an  engagement  with  m  different  types 
of  weapon  systems  on  the  X  side  and  n  for  Y  (see  Figure  11).  Although  more 
complicated  types  of  force  interactions  may  be  postulated,  we  will  consider 
the  "natural"  extension  of  (5.1.1)  to  this  combat  situation.  We  accordingly 
assume  that 

(Al)  the  attrition  effects  of  various  different  enemy  weapon-system  types 
against  a  particular  friendly  target  type  are  additive  (no  mutual  sup¬ 
port,  i.e.  no  synergistic  effects),  and 

(A2)  the  loss  rate  to  each  enemy  weapon-system  type  is  proportional  to  the 
number  of  enemy  firers  of  that  type. 

Let  Y^j  denote  those  Yj  who  engage  X.,  and  let  y^.  denote  the  corresponding 
number  of  Y^  and  similarly  for  yj.  Similar  quantities  are  analogously  de¬ 
fined  for  the  X  force.  We  observe  that  we  then  have 

(5-51) 

We  now  introduce  the  allocation  factor  =  y4i/y.  =  fraction  of  Y..  who  en- 

'  *J  'J  J  J 

gage  X^.  It  follows  that 

yij  =  *ijyj'  (5-5,Z) 

Let  a^j  denote  the  "inherent"  weapon-system  kill  rate  of  Yj  against  live  Xi 
targets,  i.e.  the  rate  at  which  one  Yj  can  kill  X.  targets. 
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X  force  (m  different  Y  force  (  n  different 

tueopon ■  system  types)  weopon- system  types) 


Figure  11:  Schematic  of  combat  between  heterogeneous  forces. 

In  this  figure  Y,j  denotes  those  Yj  who  are  engaging 
X-j,  and  y-jj  denotes  the  corresponding  number  of  Y*j 
and  similarly  for  yj.  Also,  a^j  denotes  the  "inherent" 
weapon-system  kill  rate  of  one  Yj  against  live  X, 
targets,  i.e.  the  rate  at  which  one  Yj  can  kill 
targets 


Let  us  now  examine  how  (Al)  and  (A2)  lead  to  the  following  linear  model  (no 
synergistic  effects  for  weapon  systems  in  joint  operations)  for  and  y^  >  0 
for  all  1  *  l,...,m  and  j  *  l,...,n 


dxi 

■3T 


< 


1  '*’ijaijyj  with  *i^  *  *1  f0r  '  “  1 . ra* 

j  *  1 


iE1a1jbijXi  w1th  yj*0^  =  yj  f°r  j  *  1 . n’ 


(5.5.3) 


Attrition  Modelling 


171 


Let  us  now  develop  (5.5.3)  from  assumptions  (Al)  and  (A2)  above.  Assumption 
(Al)  may  be  stated  in  mathematical  terms  as,  for  example. 


dx 

~3F 


j  n  /Xj  loss  rate\ 
f  =  '  J-1  \  due  t0  Yj/  * 


(5.5.4) 


while  assumption  (A2)  means  that 


(X.  loss  rate\  .  „ 
due  to  vj  * 


j  =  jyj ’ 


(5.5.5) 


whence  follows  (5.5.3)  from  combination  with  (5.5.4).  If  we  "absorb"  the  al¬ 
location  factors  into  the  attrition-rate  coefficients,  e.g.  let  AHi  *  ^.a.., 
then  our  linear  combat  model  (5.5.3)  may  be  written  as  (for  x.  and  y.  >  0  for 

’  J 

i  «  1 . m  and  j  =  l,...,n) 


-gr-  =  -  r  A..y,  with  x  (0)  «  x°  for  i  =  1 . m  , 

i=l  ’J  J  1  1 


dx  -j 
■ar  = 

n 

-  r 
j=l 

l< 

•nLj 

$  F 

m 

-  r 
i  =1 

(5.5.6) 


•gf  =  *  E  Bijxi  with  y.(0)  *  y.  for  j  *  1 . n 


Although  it  Is  an  easy  matter  to  develop  an  analytical  solution  to  (5.1.1) 
and  more  difficult  to  solve  (5.4.3)  with,  for  example,  the  particular  coeffi¬ 
cients  (5.4.4),  it  is  essentially  impossible  to  solve  differential  equations 
like  (5.5.6)  for  combat  interactions  with  any  degree  of  complexity.  Conse¬ 
quently,  numerical  integration  methods  are  usually  used  to  numerically  deter¬ 
mine  the  force  levels  as  a  function  of  time  x.(t)  and  y.(t)  in  complex  oper- 

1  J  g\ 

ational  models  like  BONDER/IUA.  In  such  complex  operational  models,  1  the  at¬ 
trition-rate  coefficients  A^-  and  B^,  are  (as  they  are  in  the  real  world)  com¬ 
plex  functions  of  the  weapon-system  capabilities,  target  characteristics,  dis¬ 
tribution  of  the  targets,  allocation  procedures  for  assigning  weapons  to  tar¬ 
gets,  etc.  These  models  then  attempt  to  reflect  these  complexities  by  parti¬ 
tioning  the  attrition  process  into  four  distinct  subprocesses: 

(1)  the  fire  effectiveness  of  weapon-system  types  firing  on  live  targets, 

(2)  the  allocation  process  of  assigning  weapons  to  targets. 


Essentially  all  complex  operational  Lanchester-type  combat  models  that  re¬ 
present  engagements  in  detail  (i.e.  do  not  aggregate  forces  with  firepower 
scores)  and  are  in  current  operational  use  in  the  United  States  have  been 
developed  by  the  principals  of  Vector  Research,  Inc.  The  discussion  here 
follows  that  of  Bonder  and  Farrell  (B,  pp.  11-17). 
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(3)  the  inefficiency  of  fire  when  weapon-system  types  engage  other  than  live 
targets,  and 

(4)  the  effects  of  terrain  on  limiting  firing  activities  of  weapon-system 
types  and  on  mobility  of  the  systems. 

Bonder  and  Farrell  [8,  pp.  16-17)  include  the  effects  of  the  first  three  sub¬ 
processes  above  on  the  attrition-rate  coefficient,  for  example,  as 

Aij(r)  =  ^ijIijaij(r)  •  (5'5'7) 

where  ij;.  .  denotes  the  allocation  factor  (the  fraction  of  Y,  who  are  firing  at 

y 1 J  ^ — —— — —  j 

X-),  i! .  denotes  the  intel 1 igence  factor  (the  fraction  of  Y  •  j  who  are  actual¬ 
ly  engaging  live  Xi  targets),  and  a^fr)  denotes  the  "inherent"  weapon-system 
kill  rate  (the  rate  at  which  one  Yj  kills  live  X-  targets  when  it  is  firing 
at  them).  Here,  for  simplicity,  we  have  assumed  that  the  inherent  weapon-sys¬ 
tem-kill  capability  (as  quantified  by  a^)  depends  on  only  the  range  between 
firer  and  target  (see  Bonder  and  Farrell  [8)  for  further  details).  Similar  to 
the  case  of  homogeneous  forces,  the  "inherent"  weapon-system  kill  rate  a-j  is 
computed  as 

a.  .  =  1/E  Tw  y  ,  (5.5.8) 

1J  Vj 

where  Ty-y  (a  r-v.)  denotes  the  time  for  a  single  Y,  firer  to  kill  an  X ■ 
target.  ’  J 

Thus,  Bonder  and  Farrell's  (8)  approach  (see  also  (26);  [27]  and  Cherry  [58]) 
basically  decomposes  the  battlefield  into  unit  engagements,  and  these  are  fur¬ 
ther  decomposed  into  a  series  of  one-on-one  duels  between  opposing  weapon- 
system  types.  For  each  firer-target  pair  one  must  perform  a  detailed  analysis 
of  a  single  firer  engaging  a  passive  target.  Force  interactions  are  then  tied 
together  with  attrition  equations  similar  to  (5.5.6),  and  these  assessment 
equations  are  made  to  respond  to  the  evolution  of  combat  (e.g.  changing  firer 
positions)  through  the  operational  factors  influencing  kill  rates.  Terrain  ef¬ 
fects  are  incorporated  into  such  models  by  computing  intervisibility  (i.e.  ex¬ 
istence  of  1 ine-of-sight)  for  each  target-firer  pair  based  on  their  map  loca¬ 
tions.  Consideration  is  given  to  cover,  concealment,  terrain  roughness,  etc. 
but  time  does  not  allow  us  to  go  into  further  details  here. 
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5.6  Current  Detailed  Models  of  Tactical  Engagements 

The  following  are  currently  operational  models  (used  in  the  United  States) 
that  use  detailed  lanchester-type  equations  to  assess  casualties  in  tactical 
engagements: 10 ^ 

battalion-level  combat:  BONDER/ IUA,  BONDER  A1RCAV,  BLDM,  AfISWAG,  FAST, 
division-level  combat:  DIVOPS, 
theater- level  combat:  VECTOR- 1. 

The  modern  large-scale  digital  computer  has  made  such  large-scale  models  pos¬ 
sible.  The  data  and  data-base  problems  associated  with  such  models  are,  how¬ 
ever,  formidable  although  no  less  so  than  those  for  detailed  Monte  Carlo  com¬ 
bat  simulations.  For  example,  VECTOR-1  may  require  between  200,000  and  300,000 
pieces  of  input  data  for  a  "typical"  run.  Such  models  consider  heterogeneous 
forces,  battle  plans  (ground  order  of  battle  and  air  order  of  battle),  target 
acquisition,  allocation  of  fire,  fire  support  by  ground  weapons,  movement,  in¬ 
telligence,  command  and  control,  logistics,  etc.  They  have  been  developed, 
though,  from  the  basic  analytical  structure  discussed  above  by  the  process  of 
enrichment  (which  we  have  also  considered  above). 

Although  VECTOR-1  may  well  be  the  prototype  of  the  theater-level  combat  model 
of  the  future,  as  of  August  1977  it  apparently  had  not  been  used  operational¬ 
ly  [14],  whereas  it  was  estimated  that  (for  the  same  time  period)  the  approx¬ 
imate  frequency  of  use  of,  for  example,  ATLAS  was  600  times  per  year  and  that 
of  CEM  was  25  times  per  year  (see  [I4j  for  further  details).  Even  BLDM  was 
only  used  4  to  5  times  per  year  (14).  A  partial  explanation  of  the  relatively 
low  frequency  of  use  of  such  operational  differential  combat  models  may  be 
that  they  are  rather  demanding  in  resources  (especially  highly  technically 
qualified  people  to  maintain,  modify,  and  exercise  them). 


6.  Overview  of  Aggregated-Force  Models  of  Attrition  in  Tactical  Engagements 

In  stark  contrast  to  the  detailed  Lanchester-type  models  of  attrition  in  tac¬ 
tical  engagements  are  the  aggregated-force-attrition  models  that  combine  all 
the  various  different  weapon-system  types  on  a  side  in  some  particular  geog¬ 
raphical  combat  area  (or  "sector")  into  a  single  equivalent  homogeneous  force. 


^  Documentation  on  these  models  is  discussed  in  Footnote  5);  see  also 
Bostwick  et  al.  (30),  Cordesman  [31],  and  Farrell  [32], 
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The  daily  loss  in  combat  power  as  quantified  by  the  unit's  firepower  index  is 
then  assessed  on  the  basis  of  several  operational  factors,  principal  of  which 
is  the  force  ratio  (actually  the  ratio  of  the  attacker's  firepower  index  to 
that  of  the  defender).  Current  theater-level  combat  models  typically  use 
curves  of  daily  fractional  (or  percentage)  casualties  versus  the  force  ratio 
(for  both  the  attacker  and  also  the  defender  for  each  of  several  engagement 
types  such  as  meeting  engagement,  attack  of  prepared  position,  etc.)  for  this 
assessment  of  losses.  These  curves  supposedly  have  an  empirical  basis  (see 
(1,  pp.  23-28)  or  Anderson  et  al.  [66,  p.  53);  however,  Cockrell  and  Ball 
[67,  especially  p.  1-2)  have  a  different  opinion).  Unfortunately,  there  is  no 
explicit  relationship  between  weapon-system  parameters,  operational  factors, 
and  attrition  as  there  is  for  detailed  Lanchester-type  models  (e.g.  recall 
(5.1.1),  (5.3.1),  and  (5.3.2)  above;  also  see  [261  °r  [27,  pp.  3-4)). 

Although  such  aggregated-force  models  are  much  simpler  than  the  detailed  dif¬ 
ferential  combat  models  and  therefore  more  computationally  convenient,  a 
large-scale  digital  computer  is  still  required  for  their  implementation.  Such 
aggregated-force  models  have  been  fairly  widely  criticized  (see,  for  example, 
Bonder  [7),  Honig  et  al.  [12),  or  Stockfisch  [9 )) ,  but  large-scale  conven¬ 
tional-force  ground-combat  models  that  use  such  aggregation  techniques  have 
been  and  continue  to  be  essentially  the  only  analysis  tools  used  for  large- 
scale  conventional -force  military  analyses  in  the  United  States  (see  [14)) 
and  also  NATO  countries  [33).  The  simple  fact  is  that  some  type  of  aggrega¬ 
tion  must  be  done  in  order  to  model  theater-level  combat. 


6.1  Aggregation  of  Forces  in  Combat  Analyses 

The  modern  battlefield  contains  many  diverse  weapon-system  types  that  comple¬ 
ment  each  other  and  operate  as  "combined-arms  teams".  For  example,  there  can 
be  both  mounted  and  dismounted  infantry,  tanks,  various  types  of  anti-tank 
weapon  systems,  artillery,  mortars,  infantry  with  rifles,  infantry  with  ma¬ 
chine  guns,  etc.  One  must  then  either  model  such  operations  in  great  detail 
or  find  some  means  for  aggregating  forces.  Military  planners**^  and  military 


Military  planners  have  apparently  used  the  firepower-score  approach  (see 
below  in  main  text)  for  at  least  thirty  years  (see  Mulholland  and  Specht 
[68))  to  plan  operations  and  to  plan  and  control  tactical  exercises.  Al¬ 
though  the  origins  of  using  firepower  scores  for  these  purposes  are  some¬ 
what  obscure,  they  are  still  in  use  today  (see  [69)).  Furthermore,  it  ap¬ 
pears  as  though  such  use  of  firepower  scores  in  planning  was  the  origin 
of  their  use  by  operations  researchers  in  modelling  large-scale  ground 
combat. 
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operations  analysts  have  consequently  developed  various  index-number  ap¬ 
proaches  for  aggregating  the  diverse  combat  capabilities  of  such  a  heteroge¬ 
neous  military  force  into  a  single  scalar  measure  of  combat  power.  Although 
there  are  many  such  indices*^  of  the  relative  combat  capabilities  of  mili¬ 
tary  units,  all*^  are  essentially  variations  on  the  same  theme,  and  conse¬ 
quently  we  will  generically  refer  to  any  such  index-number  approach  as  a  fire¬ 
power-score  approach. 

The  firepower- score  approach  develops  one  single  number  (referred  to  as  the 
firepower  index)  to  represent  the  "combat  potential"  of  a  military  unit.  A 
linear  model  is  used  to  develop  this  index  number,  i.e.  the  firepower  index, 
from  the  scores  of  individual  weapon  systems  as  Table  V  shows.  As  Stockfisch 
[9]  has  emphasized,  however,  the  words  score  and  index  should  not  be  regarded 
as  being  synonymous.  We  should  use  the  term  firepower  score  to  refer  to  the 
military  capability  or  value  of  a  specific  weapon  system  and  use  the  term 
firepower  index  -  which  is  obtained  by  summing  scores  -  to  refer  to  the  mili¬ 
tary  capability  or  value  of  some  aggregation  of  diverse  weapons.  In  other 
words,  the  firepower  index  of  the  X  force,  denoted  as  I x ,  is  given  by 

^  Vi  * 

where  s-  denotes  the  firepower  score  of  the  ith  X  system  and  x.  denotes  the 
number  effective  in  the  unit  (see  Table  V  again). 

Although  many  firepower-score  methods  claim  that  the  firepower  score  of  a 
weapon  system  is  determined  as  the  product  of  a  measure  of  single-round  le¬ 
thality  and  the  expected  expenditure  of  ammunition  during  a  fixed  period  of 
time,  in  actuality  varying  amounts  of  subjectivity  are  involved  in  the  devel¬ 
opment  of  such  a  firepower  score.  For  this  and  other  reasons  (e.g.  see  [12]), 
the  firepower-score  approach  has  received  a  fair  amount  of  criticism.  Never¬ 
theless,  it  is  essentially  the  only  approach  that  has  been  used  to  model 
large-scale  combat  in  currently  operational  ground-combat  models  (e.g.  see 


^  Examples  of  such  scores/indices  are  given  in  Footnote  3)  above.  Bode  [70] 
has  given  an  excellent  discussion  of  the  use  of  such  index  numbers  in  gen¬ 
eral-purpose  force  analysis,  while  Aldrich  and  Bode  [71]  have  given  a  lu¬ 
cid  discussion  of  the  conceptual  problems  of  aggregation  in  theater-level 
combat  models. 

The  one  exception  is  the  antipotential  potential  or  WEV  (see  Footnote  3) 
above,  Howes  and  Thrall  [72],  and  Anderson  [73];  [74]),  which  may  be  ex¬ 
ercised  in  the  running  of  IDAGAM  (see  Anderson  et  al .  [66]).  ATLAS  and 
other  models  that  employ  the  firepower-score  approach,  however,  are  cur¬ 
rently  much  more  widely  used  in  the  United  States  than  IDAGAM  (see  [14]). 
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Table  V:  Hypothetical  Example  of  Determination 
of  Firepower  Index  for  a  Combat  Unit 


Weapon 

dumber 

■  ■  n 

F l  repower 
Score 

_ 

Total  Contribution 
to  Firepower  Index 

Rif  1C , 

M-16,  5.56mm 

6,000 

1 

6,000 

MG, 

M-60,  . 30  cal 

150 

6 

900 

MG, 

M-2,  .50  cal 

250 

10 

2,500 

Mortar, 

M-125,  y 1  mm 

50 

20 

1,000 

tl owj  tzer , 

M-109  (SP)  ,  lSbrtur. 

50 

40 

2,000 

Uowi t  26 r , 

8H 

8 

30 

240 

Tank , 

M60A2 

200 

1AQ 

20,000 

TOTAL  FIREPOWER 

INDCX  32,640 

Firepower  Index  tor  U.  S.  Army's  ?th  Intantry  Division 


[14]).  In  other  words,  unless  one  duplicates  large-scale  combat  in  detail, 
one  must  use  some  type  of  index-number  approach  to  aggregate  the  many  differ¬ 
ent  types  of  forces  involved  in  modern  large-scale  military  operations  (see 
last  paragraph  of  Section  5.6).  Thus,  although  it  has  received  varying  amounts 
of  criticism  from  different  sources,  the  firepower-score  approach  is  used  by 
essentially  all  currently  operational  large-scale  ground-combat  models. 

In  large-scale  (i.e.  division-level  and  above)  ground-combat  models,  fire- 

141 

power  indices  are  used  as  a  surrogate  for  unit  strength  to:  ' 

(1)  determine  engagement  outcomes, 

(2)  assess  casualties,  and 

(3)  determine  FEBA  movement. 


l4)  Many  times  the  first  assessment  (i.e.  determination  of  engagement  out¬ 
come)  is  omitted.  For  example,  ATLAS  and  IDAGAfl  only  do  the  last  two  as¬ 
sessments.  However,  some  models  (e.g.  Theater  Battle  Model  (T8M-68)  (Ij) 
determine  the  outcome  of  an  engagement  (e.g.  whether  or  not  an  attack  is 
successful)  before  assessing  casualties.  In  this  case,  the  casualty-as¬ 
sessment  curves  depend  on  the  engagement's  outcome  (see  Figures  4  through 
7  of  (  1 1 ). 
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The  force  ratio  is  a  major  factor  (but  not  the  only  one)  used  to  make  such 
assessments.  Here,  however,  the  term  force  ratio  means  the  ratio  of  the  at¬ 
tacker's  firepower  index  to  that  of  the  defender.  Consider,  for  example,  the 
7th  Infantry  Division  of  the  U.S.  Army  and  assume  that  the  firepower  scores 
and  other  data  shown  in  Table  V  apply.  Then  the  7th  Infantry  Division  would 
have  a  firepower  index  of  32,640.  If  an  attacking  enemy  army  group  were  to 
have  a  firepower  index  of  146,880,  then  we  would  have  a  force  ratio  of  4.5 
(A/D),  where  A  refers  to  the  attacker  and  D  to  the  defender. 


6.2  General  Mathematical  Structure  of  Attrition  Calculations  in 
Aggregated-Force  Models 


The  usual  approach  (e.g.  see  l 10 ] )  for  assessing  casualties  in  firepower- 
score-based  combat  models  is  to  have  daily  casualties  (i.e.  the  casualty 
rates)  depend  directly  on  the  following  two  factors: 


(FI)  the  force  ratio,  and 
(F2)  the  engagement  type. 

It  will  be  instructive  for  us  to  hold  the  last  factor  constant  and  further 
examine  how  casualty  assessment  depends  on  the  firepower  scores  and  indices. 


The  basic  mathematical  structure  of  the  attrition  calculation  in  aggregated- 
force  models  may  be  thought  of  as  being  done  in  two  steps  and  may  be  explained 
as  follows: 


STEP  (I) 

(Aggregation  of  Forces) 


X  0 

x0  =  SiXi  * 


"v  Y  0 

y°  =  ial  ^  ’ 


(6  2.1) 


STEP  (II) 

(Mutual  Attrition  of 
the  Aggregated  Forces) 


(■  1  ar)  ■  #(?)  ”,tl'  *<°)  ■  "o’ 
(•  p  af)  ■ 8  (?)  “,th  J'i°>  •  "o’ 


(6.2.2) 


where  s*  denotes  the  firepower  score  of  the  ith  X  weapon-system  type,  x9  de¬ 
notes  the  initial  number  of  the  ith  X  system,  x^  denotes  the  initial  value  of 
the  firepower  index  for  the  X  force,  x(t)  denotes  its  value  at  time  t,  A(  */y) 


denotes  a  given  function  of  the  force  ratio,  t  =  0  denotes  the  start  of  the 
attrition  calculation,  and  similarly  for  the  corresponding  Y  quantities.  This 
calculation  is  then  repeated  for  each  “sector"  on  the  battlefield.  Thus,  ca¬ 
sualties  in  terms  of  a  loss  in  the  force's  combat  power  are  computed  from  some 
expression  like  (6.2.2).  In  other  words,  we  only  know  how  much  the  force's 
combat  power  was  reduced  by  a  day  of  combat  action,  and  losses  of  individual 
component  weapon-system  types  must  be  obtained  by  some  means  of  disaggrega¬ 
tion. 

ATLAS  basically  computes  casualties  in  the  above  manner,  with  the  firepower 
X  Y 

scores  (i.e.  s^  and  s^)  being  held  constant  over  time.  However,  IDAGAM  dynam¬ 
ically  recomputes  weapons'  values,  which  correspond  to  the  firepower  scores 
X  Y 

S^  and  s^  above,  according  to  the  antipotential-potential  (or  eigenvector) 
method  (see  Howes  and  Thrall  [72]  or  Anderson  [73];  [74]).  The  latter  calcu¬ 
lation  involves  the  numbers  of  enemy  targets,  allocations  of  friendly  fire, 
and  kill  probabil ities  against  enemy  targets. 

Vie  have  given  the  basic  structure  for  attrition  calculations  in  aggregated- 
f once  models  above.  In  actual  application  such  models  give  attention  to  a  mul¬ 
titude  of  details  on  combat  operations,  e.g.  positioning  of  units,  logistics 
considerations,  allocation  of  fire  (especially  supporting  fires),  air  defense, 
air  operations  including  allocation  of  aircraft  to  tactical  missions,  unit 
breakpoints,  terrain  factors,  intelligence,  command  and  control,  order  of  bat¬ 
tle,  etc.  (e.g.  see  documentation  on  CEM  [17];  [18]  or  IDAGAM  [66]  for  fur¬ 
ther  details).  Such  operational  and  tactical  factors  influence  exactly  how 
(6.2.1)  is  computed. 

6.3  Fitting  a  Differential-Equation  Model  to  Loss-Rate  Curves  typically  used 
to  Model  Large-Scale  Ground  Combat  Attrition 

In  this  section  we  will  develop  a  general  attrition  model,  whose  general  form 

fits  the  shape  of  most  loss-rate  curves  typically  used  to  model  large-scale 
151 

ground  combat.  ’  All  currently  operational  large-scale  combat  models  in  one 

Examples  of  such  casualty-rate  curves  may  be  found  in  the  documentation 
for  the  following  large-scale  ground-combat  models  (see  also  Footnote  5) 
above):  ATLAS  [10];  [16],  CEM  [17];  [18],  TBM-68  [1],  and  TAGS  [75] ;  [  76  ]. 
See  [  12 ]  for  a  general  discussion  about  such  large-scale  models  (but  for 
the  period  before  1971).  Although  IDAGAM  does  not  use  firepower  scores 
(see  Footnote  3  above),  it  uses  the  same  casual ty-rate  curves  as  ATLAS 
(see  (66,  p.  53  ).  In  fact,  it  is  stated  on  p.  53  of  [66]  that  until  bet¬ 
ter  historical  data  is  available,  the  standard  functional  relationships 
(used  in  ATLAS)  between  force  ratios  and  percent  casualties  must  still  be 
used.  Finally,  models  used  for  NATO  planning  also  employ  the  firepower- 
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way  or  another  assess  casualties  for  each  side  by  using  such  a  loss-rate  curve 
consisting  of  casualty  rate  (expressed  as  a  fraction  or  percentage  of  current 
strength  lost  per  unit  time)  plotted  against  the  force  ratio.  Here,  as  above, 
the  term  force  ratio  means  the  ratio  of  the  firepower  index  of  the  attacker 
to  that  of  the  defender,  denoted  as  A/D.  Also,  loss  here  means  loss  of  value 
for  the  side's  firepower  index,  which  can  then  be  disaggregated  into  losses 
in  numbers  of  different  weapon-system  types. 

In  other  words,  the  firepower-score  approach  takes  each  side's  heterogeneous 
forces  and  converts  them  into  an  equivalent  homogeneous  force  quantified  in 
terms  of  a  firepower  index,  daily  reduction  in  each  side's  capability  (ex¬ 
pressed  as  a  reduction  in  firepower  index)  is  then  determined  from  the  ratio 
of  the  two  such  firepower  indices,  and  finally  casualties  (i.e.  losses  in 
numbers  of  the  different  weapon-system  types)  are  assessed  by  some  means  of 
disaggregation.  We  will  now  discuss  how  a  relatively  simple  pair  of  differen¬ 
tial  equations  may  be  used  to  model  this  process  and  fit  these  loss-rate 
curves. 

Let  us  first,  however,  slightly  modify  the  equations  for  Helmbold-type  combat 
(5.1.8)  by  adding  terms  for  "operational"  losses,  i.e.  losses  not  due  to  ene¬ 
my  action  (e.g.  losses  due  to  sickness,  accidents,  etc.;  see  Taylor  and  Parry 
(49,  pp.  523-524]).  If  we  add  terms  for  such  operational  losses,  then  equa¬ 
tions  (5.1.8)  become 


1-W 

3Y  =  -a(t)-(y)  Y-y '  B(t)x  with  x(°)  =  xo  * 
1-W 

%  =  *b(t)'(x)  X'x  '  “Wy  with  y(0)  =  y0  . 


(6.3.1) 


In  equations  (6.3.1)  we  have  also  added  the  feature  of  giving  each  side  its 
own  Weiss  parameter.  In  other  words,  the  firepower-modification  factors  Ex 
and  Ev  are  no  longer  necessarily  the  same  for  both  sides,  i.e.  Ev(u;Wv)  = 

U*V,  ex(ll;Nl)  .  U‘-“X. 

For  the  case  of  constant  attrition-rate  coefficients,  (6.3.1)  becomes 


15)  cntd. 

score  approach  and  similar  casualty-rate  curves  (e.q. 
298]). 


see  (33,  pp.  287- 


-y  -  ex  with  x(0)  *  xQ  , 

(6.3.2) 

1  “6 

^  =  -b‘(j)  **  -  with  y(0)  =  yQ  , 

where  for  notational  convenience  we  have  denoted  Wy  simply  as  d  and  as  e. 
For  our  model  (6.3.2),  for  example,  X's  fractional  casualties  per  unit  time 
are  now  given  by 

/  _!  dx\  /X*s  fractional  casualties  \ 

V  x  <ftj  ~  \  per  unit  time  ) 

=  au‘d  +  e  =  avd  +  8  .  (6.3.3) 

In  Figure  12  we  show  the  relation  between  X's  fractional  casualties  per  unit 
time  and  the  force  ratio  v  =  y/x  for  the  case  in  which  X  defends  (cf.  Figure 
5).  Figure  13  shows  the  same  type  of  relation  when  X  attacks. 
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Figure  12:  Relation  between  X's  fractional  casualty 
rate  and  the  force  ratio  for  the  model 

dx  /  x  ”d 

»  -a-(-j  y  "  e*  when  X  defends 
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Figure 


13:  Relation  between  X's  fractional  casualty 
rate  and  the  force  ratio  for  the  model 

^  =  *a‘(y)  y  "  8*  when  X  attacks 


Essentially  all  of  the  principal  large-scale  ground-combat  models  currently 
in  operational  use  in  the  world  today*®)  assess  casualties  using  the  firepow¬ 
er-score  concept  and  (in  one  form  or  another)  casualty-rate  curves  of  the 
form  shown  in  Figure  14,  which  is  taken  from  documentation  on  ATLAS  [lO].Such 
casualty-rate  curves  are  typically  plots  of  fractional  casualties  per  unit 
time  (or  its  equivalent)  versus  the  force  ratio  (A/D)  for  different  engage¬ 
ment  types. Thus,  two  such  plots  like  those  shown  in  Figure  14  are  used  to 
assess  casualties,  one  curve  for  the  attacker  and  one  curve  for  the  defender. 
It  turns  out  now  that  the  Helmbold-type  model  (6.3.3)  gives  a  remarkably  good 
fit  to  almost  all  these  casualty  rate  curves,  i.e.  compare  Figures  12  and  13 
with  Figure  14  (i.e.  Figure  6-6  on  p.  6-5  of  1 10) ) ,  Figure  3  on  p.  12  of  (75], 
or  pp.  28-31  of  (76). 

In  other  words,  if  (for  a  given  engagement  type)  we  assume  that  the  fractional 
casualty  rate  depends  on  only  the  force  ratio,  then  the  so-called  (77]  asymp¬ 
totic-power  form  (6.3.3)  gives  a  very  good  fit  to  most  such  casualty-rate 


*6)  See  Footnote  14)  and  also  Footnote  5). 

*7)  For  example,  as  shown  in  Figure  14,  ATLAS  (10]  distinguishes  between  sev¬ 
en  different  types  of  engagements. 
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ATLAS 

DIVISION  CASUALTY  RATES  AS  A  FUNCTION  OF  FORCE  RATIO 
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Figure  14:  Typical  casualty-rate  curves 
used  in  ATLAS  (from  (10)) 


curves  currently  used,  and  thus  the  Helmbold-type  equations  (6.3.2)  may  be 
considered  to  model  the  attrition  process,  with  the  parameters  a,  b,  a,  e,  d, 
and  e  depending  on  the  type  of  engagement.  Moreover,  there  are  even  computer¬ 
ized  routines  available  for  the  least-squares  estimation  of  these  parameters 
(e.g.  see  [77]  especially  Figure  1  on  p.  6). 

As  we  discussed  in  Section  5.2  above,  the  model  (6.3.3),  equivalently  (6.3.2), 
can  accommodate  a  wide  variety  of  classic  attrition-rate  forms,  and  further¬ 
more  a  variety  of  attrition-rate  forms  have  indeed  been  used  in  large-scale 
ground-combat  models  over  the  years.  For  example,  ground-combat  attrition  in 
the  original  version  of  TAGS  was  assumed  to  follow  the  logarithmic  law  (see 
(78,  p.  29) ),  cf.  d  *  e  =  0  in  (6.3.2).  Today,  attrition  is  usually  modelled 
as  being  "intermediate"  between  the  logarithmic  and  square  laws.  For  example, 
comparing  Figure  12  above  to  Figure  14  (i.e.  Figure  6-6  of  [10j),  we  find  that 
the  casualty  rate  for  a  defending  force  is  best  fit  by  d  near  1  (i.e.  dx/dt  = 
-ay  -  ex).  However,  comparing  Figure  13  above  to  Figure  14  we  find  that  a 
value  for  d  around  1/2  seems  more  reasonable  for  the  attrition-rate  of  an  at¬ 
tacking  force  (i.e.  dx/dt  =  -ax^y^  -  Bx).  All  these  attrition-rate  func¬ 
tional  forms  may,  of  course,  be  handled  by  the  Helmbold-type  equations  of 
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warfare  with  operational  losses  (6.3.2)  by  taking  the  appropriate  values  for 
the  fire-effectiveness-modification  exponents  d  and  e.  Thus,  this  general 
model  (6.3.2)  has  the  flexibility  of  fitting  a  wide  variety  of  attrition-rate 
forms  that  have  been  used  to  model  large-scale  ground  combat. 

Let  us  finally  note  here  that  the  author  knows  of  no  acknowledgment  of  the 
possibility  that  the  casualty-rate  curves  such  as  we  have  been  discussing 
could  be  fit  by  a  differential -equation  model,  or  might  even  have  arisen  from 
a  formal  or  informal  understanding  of  simple  differential  equations.  Thus,  we 
have  developed  an  important  simplified  analytical  model  of  large-unit  attri¬ 
tion. 


6.4  Current  Aggregated-Force  Models  of  Large-Scale  Tactical  Engagements 

The  following  are  currently  operational  theater-level  combat  models  that  use 
the  firepower-score  approach  to  aggregate  forces  for  assessing  casualties  in 
the  manner  discussed  above: 

TAGS. 

ATLAS, 

CEM,  and 
IDAGAM. 

These  are  essentially  the  only  operational  models  currently  available  in  the 
United  States  for  analyzing  simulated  theater-level  combat.  It  was  estimated 
1 14 )  that  as  of  August  1977  the  approximate  frequency  of  use  of  ATLAS  was  600 
times  per  year,  that  of  CEM  was  25,  and  that  of  IDAGAM  II  was  between  150  and 
200. 


7.  Final  Remark! 

This  paper  has  attempted  to  survey  approaches  for  attrition  modelling  used 
today  in  the  United  States.  A  major  problem  is  that  very  little  of  this  work 
is  ever  documented  (see  Shubik  and  Brewer  [13])  (let  alone  published  in  the 
open  literature).  Consequently,  many  of  the  most  important  conceptual  model¬ 
ling  Issues  are  either  not  well  articulated  or  never  explicitly  stated  at  all 


18) 


References  to  documentation  about  these  models  may  be  found  in  Footnotes 
5)  and  15)  above. 
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The  author  believes  that  this  seminar  of  the  Hochschule  der  Bundeswehr  Miin- 
chen  is  an  important  step  in  improving  communications  among  OR  workers  inter¬ 
ested  in  combat-modelling  methodologies. 
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